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ABSTRACT 


An  improved  understanding  of  how  one  fluid  displaces  another 
in  a  porous  medium  can  be  gained  by  measuring  saturation  as  a  function 
of  time  and  distance,  a  procedure  rarely  undertaken.  To  this  end,  a 
series  of  immiscible  displacements  were  conducted,  at  various  viscosity 
ratios,  in  which  in  situ  saturations  were  measured  using  a  relatively 
new  technique  based  on  micro-wave  attenuation.  The  experimental 
saturation  profiles  obtained  in  this  manner  were  compared  with  those 
obtained  by  numerically  solving  a  recently  derived  Lagrangian 
formulation  of  the  immiscible  displacement  equation.  Finally,  the 
results  obtained  were  analyzed  in  the  light  of  published  scaling 
groups,  and  a  recently  proposed  stability  criterion.  The  results  of 
this  analysis  suggest  that  one-dimensional,  immiscible  displacement 
theory  unconditionally  represents  displacements  carried  out  at 
favourable  mobility  ratios.  For  adverse  mobility  ratio  displacements, 
in  horizontal  systems,  a  stability  criterion  must  be  satisfied  before 
the  theory  can  be  said  to  represent  the  displacement  process. 
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CHAPTER  I 


INTRODUCTION 


1.1  BACKGROUND 

The  macroscopic  law  governing  flow  of  a  liquid  in  porous 
media  was  developed  by  Darcy  in  1856.  In  horizontal,  one-dimensional 
systems,  Darcy's  law  can  be  written  as: 


K  dP_ 
u  dx 


(1-1) 


The  macroscopic  measure  of  resistance  to  flow  for  a 
particular  porous  medium  is  the  absolute  permeability,  K,  which  must 
be  determined  experimental ly.  Absolute  permeability  of  a  porous 
medium  is  invariant  with  the  physical  properties  of  the  liquid  used 
to  evaluate  it. 

In  all  hydrocarbon  bearing  porous  media,  it  is  well  known 

that  water  co-exists  with  hydrocarbons  which  are  either  in  the  liquid 

or  vapour  phase.  When  water  flows  with,  or  displaces,  these  hydrocarbons , 

Darcy's  law,  as  formulated  for  the  flow  of  a  single  fluid,  is  not 

applicable.  The  concept  of  absolute  permeability  is  still  valid 

because  it  is  a  property  of  the  porous  medium.  To  study  the  flow  of 

55 

hydrocarbons  with  water  through  porous  media,  Wyckoff  and  Botset 
[1936]  used  gas-water  and  oil-water  mixtures  to  flow  through  an 
unconsolidated  sand  pack.  Their  results  indicated  that  the  flow  of 
either  fluid  can  still  be  predicted  by  a  Darcy-type  law,  provided 
the  absolute  permeability  is  replaced  with  an  effective  permeability 
to  that  fluid.  Effective  permeability  is  a  property  of  a  particular 


1 


. 


2 


rock  and  is  also  a  function  of  the  relative  amount  of  the  fluid 
present  in  an  elemental  volume  of  the  porous  medium.  This  relative 
amount  is  the  fraction  of  pore  volume  occupied  by  a  fluid  in  an 
element  of  the  porous  medium  and  is  known  as  the  saturation  of  that 
fluid.  Darcy's  law  for  water  flowing  together  with  hydrocarbons  is 


vw  =  - 


WV  3Pw 


(1-2) 


Effective  permeabilities,  although  unique  functions  of 
saturation  for  a  fluid  and  a  porous  medium,  do  have  implicit  in  them 
effects  of  interfacial  tensions  between  the  fluids  and  the  wettability 
characteristics  of  the  porous  medium  for  those  fluids.  They  also 
depend  on  saturation  history. 

The  current  immiscible  displacement  theory  is  based  on 
Darcy's  law  { Eq . ( I -2 ) }  written  for  both  the  fluids,  an  equation  to 
link  the  two  flow  relationships  and  conservation  of  volume  in  the 
porous  medium.  The  linkage  between  the  two  flow  relationships  is  the 
defining  equation  for  capillary  pressure  [Leverett  ,  1941].  This 
formulation  has  inherent  uncertaini ties  as  it  requires  experimentally 
determined  effective  permeabi 1 ity-saturation  and  capillary  pressure- 
saturation  curves.  More  often  than  not  mathematical  models  are  used, 
rather  than  experimental ly  derived  relationships. 


1.2  PURPOSE  OF  THE  STUDY 

The  conventional  immiscible  displacement  theory  gives  rise 
to  a  second  order,  non-linear,  partial  differential  equation  [Rapoport 
and  Leas^,  1953]  which  has  no  known  analytical  solution.  Despite 
uncertainties  in  the  relative  permeability  and  capillary  pressure 
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curves,  the  immiscible  displacement  equation  is  solved  numerically, 
resulting  in  one-dimensional  saturation  profiles.  It  was  the  purpose 
of  this  study  to  experimentally  measure  these  saturation  profiles  and 
then  analyse  them,  quantitatively,  together  with  the  numerical 
solutions.  These  numerical  solutions  were  generated  using  published 
relative  permeability  and  capillary  pressure  models. 
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Rose  [1974,  p . 233]  recognized  the  dearth  of  this  type  of 
analysis  in  the  literature,  and  mentions  that 

"in  the  last  decade,  much  progress  has  been  made  in 
the  development  of  stable  and  rapidly  convergent 
numerical  solution  techniques  for  the  second  order 
non-linear  partial  differential  equation.  This  work 
has  yielded  at  least  qualitatively  plausible  forecasts 
of  the  characteristics  of  immiscible  displacement 
processes.  Still  no  definitive  laboratory  work  has 
been  reported  in  the  literature  by  which  the 
quantitative  significance  of  the  theory  can  be  judged". 

Moreover,  he  also  suggested  that,  to  ascertain  the  consistency  of  the 

theory,  saturations  should  be  measured  as  a  function  of  time  and 

distance  as  this  approach  was  generally  considered  superior  to  the 

analysis  of  the  total  pressure  drop  versus  time  curve,  which  has  been 

53 

carried  out  by  Warren  and  Calhoun  [1954]. 

1.3  METHODOLOGY 

The  following  methodology  was  devised  and  executed. 

Saturations  were  measured  as  a  function  of  time  and 
distance  employing  a  suitable  technique,  selected  on  its  merits,  as 
enumerated  in  the  literature.  Displacement  tests  were  conducted  at 
different  viscosity  ratios  in  an  unconsolidated  sand  pack.  Rate  was 
also  a  parameter.  A  suitable  published  numerical  procedure  was 
chosen  to  generate  calculated  saturation  profiles.  A  numerical 
procedure  was  also  developed  in  this  study.  The  procedure  which  gave 


, 
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a  closer  match  to  the  experimentally  measured  saturation  profiles  was 
used  for  comparison  purposes.  For  support  of  numerical  procedures, 
capillary  pressure  and  relative  permeability  models  were  taken  from 
the  literature. 

Finally  an  overall  analysis  of  the  experimental  and  numerical 
results  was  envisaged  in  the  light  of  established  scaling  groups  and 
stability  criteria. 


CHAPTER  II 


THEORY  AND  LITERATURE  REVIEW 

11. 1  IMMISCIBLE  DISPLACEMENT 

11. 1 .1  Equations 

Immiscible  displacement  theory  is  used  to  predict  the 

saturation  behaviour  during  displacement  of  oil  by  water.  This 

theory  is  based  on  the  solution  of  the  following  equations: 
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Darcy's  law  as  written  for  both  oil  and  water  , 
[Wyckof  and  Botset,  1936] 
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the  capillary  pressure  definition  ,  [Leverett,  1941] 

P  (S  )  =  P  (S  )  -  P  (S  ) 

cv  wy  ov  w'  wv  viJ 


( 1 1-3) 


and  the  continuity  or  material  balance  equation  for  incompressible 
fluids^0,  [Buckley  and  Leverett,  1942] 
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Rapoport  and  Leas47  [1953]  combined  these  four  equations 
to  illustrate  that,  according  to  this  theory,  the  physics  of 
immiscible  displacement  in  porous  media  can  be  described  by  the 
following  second  order,  non-linear,  parabolic,  partial  differential 
equation 
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Bentsen  [1976]  nondimensional ized  these  equations  [  ( 1 1 - 1 ) 
(II-4)]  utilizing  the  procedure  of  Heliums  and  Churchill27  [1961]. 

The  normalized  equations  when  written  for  horizontal  systems  are, 
af. 
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G ( S* )  =  Ftl(S*)  (for  horizontal  systems) 
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Normalization  of  capillary  pressure,  P  ,  and  effective 
permeabilities,  KeQ  and  K0W,  required  special  treatment. 
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II. 1.2  Normalization  of  Capillary  Pressure  and  Effective 
Permeabilities 

Bentsen  and  Anli^’3  [1976,  1977]  proposed  a  capillary 
pressure  model  and  used  data  obtained  from  a  centrifuge  to  validate 
it,  as  well  as  to  estimate  the  parameters  postulated  in  the  model. 
This  model  is  for  a  drainage  capillary  pressure  versus  saturation 
curve  but  for  the  purpose  of  this  study  it  is  also  considered  to 
represent  an  imbibition  capillary  pressure  versus  saturation  curve. 
The  model  is , 


or 
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S* 


(11-16) 


(II-17) 


where 

Pr 

=  -£•  (Il-li 

C  q 

and  the  capillary  normalizing  parameter,  a,  has  the  units  of 

pressure.  Mathematically,  a  is  the  area  under  the  tt  versus  S*  curve 

and  physically  it  is  influenced  by  interfacial  tension,  pore-size 

distribution  and  wettability.  Figure  1 1 - 1  shows  the  Pc  versus  Sw 

and  tt  versus  S*  curves, 
c 

Bentsen3  [1976]  normalized  the  effective  permeabilities 
after  the  manner  of  Perkins  and  Collins^  [I960].  That  is. 
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Normalized  Capillary  Pressure  (tt  )  Capillary  Pressure  (P  ) 
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CAPILLARY  PRESSURE-SATURATION  RELATIONSHIPS 

(Bentsen,  1976) 


FIGURE  II-l 
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where  K  and  K  are  effective  permeabilities  to  water  at  (1  -  S  ) 
wr  or  r  v  or 

and  to  oil  at  S  . ,  respectively.  Figure  I 1-2  shows  how  the 

W I 

normalization  of  effective  permeabilities  is  carried  out. 

II. 1.3  Review  of  Theory  and  Numerical  Solutions 

Terwilliger  et  al.^  [1951],  Welge^4  [1952]  and 

Rapoport  and  Leas^  [1953], with  their  experimental  and  analytical 

treatment  of  the  immiscible  flow  equations, enhanced  the  physical 
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understanding  of  such  displacement.  Welge  ,  essentially  using  the 

Buckley  and  Leverett^  [1942]  equations,  demonstrated  that  balancing 

areas  on  a  S*  versus  £  plot,  to  get  around  the  absurdity  of  multiple 

saturations  at  one  point,  is  equivalent  to  drawing  a  tangent  to  the 

F  (S*)  versus  S*  curve  from  S*  =  0.  [Figure  I 1-3] . 
w 

As  a  consequence  of  the  philosophical  argument  put  forward 

by  Cardwell  ^  [1959],  regarding  the  multi-valued  solutions  obtained 

using  non-capillary  Buckley-Leverett  theory,  doubts  were  cast  on  the 

experimental ly  determined  F, ,(S*)  versus  S*  curve.  Because  F  (S*) 

needs  experimentally  determined  effective  permeabilities,  Cardwell 

inferred  that  either  the  effective  permeabilities  are  erroneously 

measured  or  the  theory  is  incomplete.  During  all  this  time 

Buckley-Leverett^  theory  was  being  put  to  use  in  solving  field 

waterflooding  problems.  [Dardaganian^,  (1958);  Kufus^,  (1959)]. 

1 6 

With  the  advent  of  digital  computers,  Douglas  et  al . 
[1958]  solved  the  non-linear  partial  differential  equation  [Eq .  ( 1 1-5 ) ] 
by  a  transformation  of  the  dependent  variable,  Sw,  and  by  employing 
a  set  of  boundary  and  initial  conditions  which  allowed  for  a  mobile 
water  saturation  ahead  of  the  front.  Outlet  end  effects  were 
accounted  for  by  allowing  the  outlet  end  saturation  to  build  up  to 
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Effective  Permeability 


STEPS  IN  NORMALIZING  EFFECTIVE  PERMEABILITY 
(Perkins,  et.  al.,  1960;  Bentsen,  1976) 


STEP  1 


1 


STEP  3 


Normalized  Saturation 
(S*) 


FIGURE  1 1-2 


Normalized  Saturation 
(S*) 


n 


Buck! ey-Leverett  Solution  (1942) 


Welge's  Tangent  Construction  (1952) 


FIGURE  1 1-3 
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its  maximum  value,  1  -  SQr,  prior  to  water  breakthrough. 

Interest  was  also  focussed  on  solving  Eq.  ( 1 1 -6 )  and 

Eq.  (II-7),  either  combined  together,  [Fayers  and  Sheldon^,  1959] 

38 

or  solved  seperately,  [McEwen  ,  1959]  in  the  fractional  flow 

domain,  viz.  f  (S*,  t*).  These  studies  demonstrated  that  when  the 

w 

complete  immiscible  displacement  equation  is  solved  (i.e.  including 
the  capillary  term)  the  multi-valued  saturations  were  eliminated. 

Hovanessian  [1961]  solved  the  immiscible  displacement 

1  c 

equation  using  a  transformation  similar  to  that  of  Douglas  et  al. 
and  found  that  the  magnitude  of  the  initial  water  saturation 
distribution,  essential  in  getting  the  numerical  procedure  started, 
did  not  effect  the  shape  of  the  profiles  later  on  in  the  life  of  the 
flood. 

After  these  initial  numerical  solutions,  stress  was  placed 
on  optimizing  the  various  algorithms  for  accuracy  and  computer 
processing  time.  [Gottfried33  et  al.  (1966);  Blair3  et  al .  (1969)]. 

Numerical  solution  of  the  immiscible  displacement  equation 
in  the  saturation  domain,  S*U,  t*),  without  any  transformation,  was 
carried  out  by  Flock^  et  al.  [1977]. 

Bentsen3  [1977]  solved  Eq.  ( 1 1 -6 )  and  Eq.  ( 1 1 - 7 )  in  the 
fractional  flow  domain,  ftl(S*,  x*),  with  a  more  appropriate  inlet 
boundary  condition  which  allows  the  saturation  to  build-up  to  a 
maximum,  in  contrast  to  all  the  previous  solutions  which  fixed  the 
saturation  at  the  maximum.  This  was  physically  realistic. 

[Perkins43,  1957].  It  was  also  demonstrated  that  the  Buck! ey-Leverett1 0 
solution  is,  in  fact,  the  steady  state  solution  to  the  immiscible 
displacement  equation.  Eq . ( I 1-6)  and  Eq. (II-7),  when  combined  and 
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the  variable  £  eliminated,  results  in,  [Bentsen  ,  p.26]s 

2  3% (S*,  T*)  3 f  ( s* ,  T*) 

[f  (S*,  T*)  -  G(S*)]  - * -  =  N  C(S*)  -  -  — y  ■« -  (11-21 ) 

w  35*2  C  3t 

C 

As  a  consequence  of  Eq. (11-21),  Bentsen  postulated  that 

the  f  (S*,  x*)  function  approaches  the  G(S*)  function  in  the 
w 

fractional  flow  domain.  [Fig. I 1-4] -  The  G(S*)  function  is  the 
non-capillary  fractional  flow  curve.  At  that  instant  when  f  (S*,  T*) 
becomes  invariant  with  time  (steady  state),  Eq. (11-21)  can  be  written 
as, 

2  32f 

[ f w ( S* )  -  G(S*)]  ^  =  0  (11-22) 

For  Eq. (11-22)  to  hold,  either  of  the  terms  on  the  left 
hand  side  can  be  equated  to  zero  resulting  in  either 

f  (S*)  =  G  ( S* )  S*  <  S*  <  1  (11-23) 

or 

f  (S*)  =  f  .  U  0  <  S*  <  S*  (11-24) 

wv  ‘  wf  S|  f 

Eq. (11-23)  and  Eq. (11-24)  are  also  evident  from  Fig. 1 1 - 4 . 

g 

As  a  consequence,  the  analysis  undertaken  by  Bentsen 

in  the  f  (S*,  x*)  domain  resulted  in  an  interesting  observation  in 
w 

the  S*U,  x*)  domain.  As  the  f  (S*,  x*)  function  approaches  G(S*), 
the  saturation  profiles  tend  toward  a  fixed  configuration  for 
saturations  less  than  S$,  which  is  the  "stabilized  zone"  postulated 
by  Terwilliger51  et  al.  [1951]  and  further  analysed  by  Jones-Parra 
and  Calhoun^  [1953]. 
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BEHAVIOUR  OF  f,,(S*,  t*) 

w 

(after  Bentsen,  1977) 


FIGURE  1 1-4 
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II. 2  ROLE  OF  SATURATION  PROFILES 

It  is  well  established  that  the  solutions  to  Eq . ( I 1-6)  and 
Eq . ( I 1-7 )  or  Eq. ( I 1-5)  are  saturation  profiles  which  in  turn  dictate 
the  oil  recovery  from  porous  media.  To  solve  these  equations, 
independently  determined  or  modelled  relative  permeability  and 
capillary  pressure  functions  are  required.  The  numerical  solution 
of  the  immiscible  displacement  equation  can  be  analysed  for  different 
displacement  rates  or  viscosity  ratios.  Parallel  to  the  numerical 
studies,  displacement  experiments  can  also  be  conducted.  Various  oils 
can  be  used  to  study  the  effect  of  interfacial  tension.  By  changing 
the  type  of  porous  medium  the  role  of  wettability  can  also  be  analyzed. 

A  rigorous  test  of  immiscible  displacement  theory  would  be 
to  experimentally  measure  the  saturations  as  a  function  of  time  and 
distance.  These  saturation  profiles  could  then  be  studied  for  different 
displacement  tests. 

Two  approaches  for  analyzing  saturation  profiles  are: 
i)  To  compare  the  experimental  saturation  profiles  amongst 
themselves  by  varying  any  one  of  the  parameters. 

(viz.  rate). 

ii)  To  compare  the  experimental  saturation  profiles  with 
the  numerically  generated  ones  based  on  the  current 
immiscible  displacement  theory. 

The  second  approach  requires  a  judicious  selection  of  relative 
permeability  and  capillary  pressure  models. 
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1 1. 3  SATURATION  MEASUREMENT  TECHNIQUES 


Saturations  have  been  measured  over  a  gross  elemental 
volume,  using  the  principle  of  discrimination  between  resistivities 

cc 

of  brine  and  hydrocarbons.  Wyckoff  and  Botset  [1936]  and 
Leverett^  [1939]  used  this  technique  to  determine  relative  permeability 
curves.  Levine  [1954]  published  one  set  of  saturation  prifiles  using 
the  same  technique.  This  particular  displacement  test  was  an  oil -flood 

r  n 

at  a  favourable  viscosity  ratio  (y  /y  =  0.52).  Steber  and  Marsden 

W  0 

[1956]  suggested  one  based  on  matched  refractive  indices.  This  was  a 
step  towards  optical  techniques. 

Again  using  the  resistivity  technique,  Perkins^  [1957] 
reported  two  sets  of  saturation  profiles,  one  at  a  low  rate  and 
another  at  a  high  rate.  These  two  runs  were  waterfloods  in  the 
presence  of  connate  water.  The  viscosity  ratio  was  slightly  adverse 

(yo/pw  =  2-0)- 


Bail  and  Marsden''  [1957],  using  a  resitivity-brine  saturation 
correlation,  reported  six  runs  with  rates  ranging  from  0.00145  cc/sec 
(5.2  cc/hr)  to  0.0098  cc/sec  (35.28  cc/hr)  giving  a  superficial 

r  _  r 

velocity  for  their  core  from  7.727  (10  )  cm/sec  to  52.2  (10  )  cm/sec. 

In  their  water-floods  they  maintained  a  constant  viscosity  ratio 

(y  /y  )  of  2.44  at  30°C. 
v  0  w 

Scrutinising  the  saturation  profiles  published  upto  this  time 


it  was  evident  that  much  more  experimentation  needed  to  be  done  both 
qualitatively  and  quantitatively.  Efforts  were  being  made  in  this 
direction  and  different  techniques  were  proposed. 

X-rays  can  be  used  in  the  determination  of  liquid  saturations 
as  their  intensity  loss  exhibits  a  linear  relationship  with  saturation 
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on 

on  a  semi-log  plot.  [Geffen  and  Gladfelter  ,  1952],  For  X-rays  to 
be  effective,  the  water  has  to  be  doped  with  sodium  iodide  or  iodo- 
benzene  has  to  be  added  to  the  liquid  hydrocarbons. 

In  his  paper  on  displacement  stability,  Hagoort  [1974] 
presented  one  profile  from  an  experiment  conducted  at  a  low  rate  and 
compared  it  with  the  non-capillary  Buckley-Leverett10  solution.  The 
technique  he  used  was  based  on  selective  absorption  of  light.  He 
used  a  transparent  model  and  added  1%  (wt.)  cobalt-chloride  to  the 
water  to  colour  it  blue.  The  viscosity  ratio  in  his  experiments  was 
21.58.  He  did  not  present  any  profiles  for  those  experiments  in  which 

displacements  were  termed  "unstable". 

fn 

[1975],  taking  advantage  of  selective  absorption 

of  electro-magnetic  radiation  at  certain  frequencies  for  some 

substances,  chose  microwaves  at  21.2  GHz  which  exhibit  a  strong 

absorption  preference  to  water.  He  presented  saturation  profiles 

which  were  superior  in  quality  to  any  published  before.  He  presented 

two  sets  of  profiles,  an  oil-flood  with  a  viscosity  ratio  (y  /u  )  of 

(1/8.5)  [Parsons  ,  1978]  and  a  subsequent  water-flood  on  the  same 

system  after  connate  water  was  established.  Using  this  technique  to 

generate  saturation  profiles  during  a  micellar  flood,  Parson  and 

Jones41  [1977]  developed  the  concept  of  linear  scaling.  Baker2  [1975] 

reported  his  results  on  the  use  of  microwaves  during  a  steam-flood. 

He  reported  difficulty  in  identifying  water  vapour  in  this  kind  of 

flood  but  suggested  that  quantitative  calibration  was  possible  for 

21 

air-water,  oil-water  and  oil-alcohol  systems.  Gladfelter  [1978] 
also  used  a  set  up,  similar  to  that  of  Parsons40,  for  saturation 


measurements. 
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1 1. 4  CHOICE  OF  A  SATURATION  MEASUREMENT  TECHNIQUE 

To  generate  saturation  profiles  for  this  study,  a  technique 
based  on  microwave  attenuation  was  chosen.  Reasons  for  the  selection 
of  this  technique  were: 

i)  Resistivity  techniques  are  cumbersome  and  measurements 
are  made  at  discrete  points  with  electrodes  in  contact 
with  the  core  holder. 

ii)  X-rays  and  optical  techniques  need  doping  of  one  of  the 
liquids  with  a  foreign  substance, 

iii)  Dissection  or  post-mortem  of  the  porous  medium 

essentially  gives  a  qualitative  view  at  one  point  in 
time. 

iv)  The  microwave  technique  does  not  require  physical 

contact  with  the  core-holder  and  allows  continuous 

40 

monitoring  of  saturations  . 

v)  Among  the  published  saturation  profiles,  those 
measured  using  microwave  techniques  were  of  high 
qua! ity. 


CHAPTER  III 


MICROWAVE  THEORY,  EXPERIMENTAL  EQUIPMENT 

AND  PROCEDURE 


III.]  MICROWAVE  THEORY 

The  broad  spectrum  of  electro-magnetic  radiation  encompasses 

frequencies  ranging  from  X-rays  having  frequencies  of  more  than 

3(10^)  Hz  to  very  low  frequencies.  The  electro-magnetic  radiation 

termed  microwaves  has  frequencies  ranging  from  3(1 0^)  Hz  to  3(10^)  Hz. 

Visible  light  has  a  frequency  range  of  0.43(10^)  Hz  to  1.33(1016)  Hz, 

qa 

greater  than  microwaves  but  less  than  X-rays.  Vankoughnett 

[1973,  p . 21 ]  contends  that 

"the  fact  that  visible  light  and  microwaves  are 
both  electro-magnetic  waves  suggests  that  they 
share  common  properties." 

The  laws  governing  microwave-material  interaction, 

specifically  microwave-water  interaction,  are  important  in 

understanding  the  experiments  conducted  to  generate  saturation 

profiles.  Microwaves  interact  differently  with  different  types  of 

substances.  Basically  there  are  three  classes  of  material: 

[Tinga  ,  1970]  conductors,  insulators  and  dielectrics.  Conductors 

reflect  microwaves  so  they  are  used  to  contain  and  guide  them. 
Insulators  reflect  and  transmit  microwaves  without  any  appreciable 

loss  of  power.  They  are  almost  "transparent"  to  microwaves. 

Dielectrics,  like  water,  are  substances  which  absorb 

microwaves.  This  absorption  of  microwaves  can  be  quantitatively 

predicted  by  two  different  approaches. 
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III. 1.1  Electro-magnetic  Approach 

In  a  non-magnetic  material,  the  absorption  of  microwaves 
depends  on  the  dielectric  constant  for  that  material.  A  high 
dielectric  constant  indicates  a  greater  absorption  of  microwaves. 

A  water  molecule  has  a  dipole  moment  associated  with  it. 

When  an  oscillating  field  is  applied  to  water  molecules  they  try  to 
align  themselves  with  the  continuously  changing  field.  This 
phenomenon  is  known  as  polarization.  As  a  consequence  the  potential 
of  the  electric  field  is  reduced.  This  reduction  is  dictated  by  the 
dielectric  constant,  e',  for  the  material.  The  dielectric  constant 
is  determined  experimentally.  [Hasted  ,  et  al . ,  1953]. 

In  addition  to  the  decrease  in  potential,  there  is  a  loss 
of  power  or  irreversible  degradation  of  energy  manifested  as  heat. 

This  loss  of  power  or  "attenuation"  depends  on  the  magnitude  of 

pc 

another  experimental ly  determined  quantity,  e",  known  as  the 
"loss  factor". 

A  complex  dielectric  constant  then  can  be  defined  as, 

e*  =  e1  -  ie"  (HI-1) 

where  variation  of  e'  and  e"  with  frequency  is  shown  in  Fig.III-1. 

The  loss  factor,  e",  shows  a  definite  peak  at  a  specific  frequency 

(Fig.  111-2} .  For  water  molecules  this  peak  occurs  at  about 
52 

20  GHz  . 
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LOSS  FACTOR  VERSUS  FREQUENCY  FOR  DIFFERENT  SUBSTANCES 


FIGURE  1 1 1-2  (ref:  Parsons,  p.304) 
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III. 1.2  Optical  Approach 

Both  the  electro-magnetic  and  the  optical  treatment  of 
microwave-dielectric  interaction  is  necessary  to  the  understanding 
of  the  measurement  of  water  saturations.  Two  laws  of  radiation 
absorption,  extensively  used  in  spectroscopy,  are  quoted  here. 
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Lambert's  Law  [Harrison  ,  p.363,  1948]  which  states  that 

When  a  parallel  beam  (plane  wave)  of  monochromatic  radiation 
passes  recti  1 i nearly  through  a  homogeneous  medium,  the 
power  intensity  of  the  beam  is  reduced  as  follows, 

Ix  =  IQ  exp  { -Kx)  ( 1 1 1-2) 

where 


I  =  incident  power 

I  =  power  after  the  beam  has  traversed  a  distance,  x, 

A 

and 

K  =  extinction  coefficient. 

Beer's  Law  [Harrison  ,  p.364,  1948]  which  states  that 

"If  the  absorbing  medium  is  a  substance  in  solution,  the 
attenuation  of  radiation  on  traversing  a  given  path  length 
depends  on  the  concentration  of  the  solution". 

Beer's  Law  can  be  written  as, 

K  =  kc  (II 1-3 ) 

where 


k  =  extinction  coefficient  per  unit  of  concentration, 
and 

c  =  concentration. 

It  is  worthwhile  to  mention  here  that  the  concentration, 
c,  is  analogous  to  the  water  saturation  in  that  elemental  volume 
which  is  being  traversed  by  a  plane,  monochromatic  beam  of  radiation. 
The  remaining  substances,  like  oils  used  in  this  study,  Ottawa  sand 
(silica)  and  PVC  (material  of  core-holder) ,  are  almost  transparent 
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to  microwave  radiation  at  27  GHz.  Figure  III-2,  taken  from 
Parsons40  [p.304],  indicates  the  relative  magnitudes  of  loss  factors, 
e",  for  water  and  some  organic  compounds. 

Equations  ( 1 1 1 - 2 )  and  ( 1 1 1 - 3 )  result  in  the  relationship 
known  as  Beer-Lambert 1 s  Law, 

Ix  =  I0  exp  {-kcx}  ( 1 1 1-4) 

Several  variations  of  this  equation  are  employed,  depending 
on  the  context  of  its  usage.  For  use  in  the  discussion  of  absorption, 
as  derived  from  classical  electro-magnetic  theory  [Harrison  ,  p.369], 
the  following  equation  is  commonly  used. 

Ix  =  IQ  exp  { -  (4-rrfxK) /c}  ( 1 1 1-5) 

where 


f  =  frequency,  Hz, 
x  =  distance  traversed,  cm, 


c  =  velocity  of  light,  cm/sec. 

The  extinction  coefficient,  K,  is  related  to  the  unit  extinction 
coefficient  and  concentration  as  indicated  by  Eq . ( 1 1 1 -3 ) .  In  the 
context  of  the  experiments  conducted  in  this  study,  the  concentration 
of  the  absorptive  substance  in  the  path  of  the  microwave  beam  is 

c  =  4>S  ( 1 1 1 -6 ) 

w 

Also  x  is  the  thickness  of  the  core,  h.  One  may  substitute  for  K 
in  Eq.  from  Eq.(III-3)  and  for  c  from  Ea .  (1 1 1-6)  to  obtain 

I  =  IQ  exp  { - (47rf hkcfSw)/c}  (III-7) 


or 


I  47rfhkcj)Si 

,  x  _  _  w 

log10  I  “  2.303c 
o 


(II 1-8 ) 


, 
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III. 1.3  Relationship  Between  Electro-magnetic  and  Optical  Approaches 
The  two  approaches  can  be  related  and  the  analysis  by 
Dekker^5  [1957,  p . 1 56]  is  used  to  illustrate  this.  When  there  is 
absorption,  the  electrical  component  of  the  "light  wave"  polarized 
in  the  y-direction  and  propagated  in  the  x-di recti  on  may  be 
represented  as 

E  (x,  t)  =  E.  exp  {-(27rfkx/c)}  exp  {i2rf(t  -  —■)}  (HI-9) 

VI  ^ 

where 

E.  =  electric  field  strength  at  x  =  0, 

E  (x,t)  =  electric  field  strength  at  a  distance  x, 

y 

after  time  t, 


and 

n  =  refractive  index. 

This  equation  can  be  written  as 

E  (x,t)  =  E.  exp  { i2irf (t  -  n*x/c)}  (III-10) 

where  the  complex  index  of  refraction  is 

n*  =  n  -  i k ,  (1 1 1-11) 


For  a  non-magnetic  insulator,  the  velocity  of  propagation 
is  given  by 

v,  =  c//r~,  (in-12) 

<P 

while  the  index  of  refraction  is  defined  as 

n  =  c/v  (HI-13) 

Equations  (III-12)  and  (III-13)  lead  to 

e.  =  n2  (HI-14) 

which  can  be  written  for  a  dielectric  as 

e*  =  n*2 


(III-15) 


. 
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Re-writing  Eq.(III-l)  and  Eq.(III-n), 

£*=£*-  1e"  (III-l) 

and 

n*  =  n  -  ik  (III-ll) 

Substituting  for  e*  and  n*  in  Eq. (III-15)  yields 

e 1  -  ie"  =  n2  -  2ink  -  k2.  (III-16) 

Equating  the  real  and  imaginary  parts  one  obtains 

e'  =  n2  -  k2  (III-17) 

and 

e"  =  2nk  (III-18) 

These  equations  relate  the  extinction  coefficient,  k,  to  dielectric 
constant,  e‘,  and  the  loss  factor,  e". 


II 1 . 1 . 4  Quantitative  Evaluation  of  the  Extinction  Coefficient 

For  the  analysis  of  microwave  responses  obtained  in  this 
work,  a  theoretical  guideline  is  given  by  Eq. (III-8).  This  equation 
indicates  that  for  a  constant  porosi ty-thickness  product,  <j>h,  along 
the  length  of  the  core,  a  plot  of  log  I  versus  Sw  will  result  in  a 
straight  line  at  a  fixed  frequency,  f.  In  order  to  do  a  preliminary 
analysis  so  as  to  get  some  feel  for  this  correlation,  the  equation 
can  be  solved  for  varying  power  inputs,  I  ,  and  water  saturations, 

S  ,  to  have  an  estimate  of  attenuated  power  outputs,  I.  To  this  end 
firstly,  a  theoretical  estimate  of  the  extinction  coefficient,  k,  is 
needed. 

Hasted,  et  al .  have  determined  the  dielectric  constant, 
e',  and  the  loss  factor,  e",  experimentally,  for  water  at  20°C  and 
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a  microwave  frequency  of  23.6  GHz.  Using  Eq . ( 1 1 1-1 7)  and  Eq.(III-18) 
and  the  experimentally  determined  values  of 
e1  =  30.5 

and 


e"  =  35.6 

a  value  for  the  extinction  coefficient,  k,  is  obtained  as  2.86.  A 
simple  analysis,  using  this  theoretical  value  of  k,  is  undertaken  in 
Appendix  C. 

III. 2  MICROWAVE  INSTRUMENTATION 
III. 2.1  Antennas 

Special  kinds  of  directional  antennas  are  needed  to  have  a 
linearly  polarized,  plane,  wave  of  electro-magnetic  radiation.  The 
transmitting  and  receiving  antennas  have  to  be  identical  by  the 
general  theorem  of  electro-magnetic  reciprocity.  [Chu  and  Barrow  , 
1939,  p.333] . 

For  the  purpose  of  measuring  water  saturations  a  pyramidal 

40 

horn  with  a  lens  antenna  to  alter  the  radiation  pattern  was  selected  . 
The  pyramidal  horn  transmits  a  spherical  wave  with  primary  and 
secondary  lobes,  but  when  a  lens  is  placed  on  the  mouth  of  the  horn 
this  spherical  wave  is  transformed  to  a  plane  wave,  given  a  suitable 
choice  of  refractive  index  and  shapes  of  lens  surfaces.  [Fradin  , 
1961,  p. 13] . 

Based  on  the  theory  of  refraction,  in  optics,  [Fradin  , 
p . 265]  a  pair  of  lenses  were  designed  for  a  pyramidal  horn.  The 
details  are  given  in  Appendix  C.  Polystyrene  was  chosen  as  the 
material  for  the  lens  because  of  its  low  loss  factor,  e".  The 


. 
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refractive  index  of  the  polystyrene  used  was  1.6  [Brown9,  1953,  p.32]. 

111. 2. 2  Instrumentation 

A  schematic  plan  view  of  the  experimental  set-up  is  shown 
in  Figure  III-6.  Microwaves,  at  a  constant  frequency  of  27  GHz, 
were  generated  by  an  osci llator/R. F.  unit.  (Hewlett  Packard,  HP8679A). 
This  unit  has  an  automatic  level  control  which  senses  a  feedback  from 
a  crystal  detector  (HP  R422A).  The  crystal  detector  takes  a  minute 
fraction  of  the  radiation  with  the  help  of  a  directional  coupler 
(HP  R752D)  placed  downstream  of  the  oscillator.  A  filter  (HP  R362A) 
and  a  tuner  were  also  placed  in  the  circuit  between  the  directional 
coupler  and  the  oscillator.  The  operating  frequency  of  27  GHz  lies 
in  the  R-band  of  the  microwave  spectrum,  so  appropriate  waveguides 
were  used  for  transmission.  From  the  directional  coupler  the 
microwaves  were  transmitted  to  the  sending  horn  and  lens  antenna. 

The  sending  and  receiving  horns  are  shown  in  position,  relative  to 
the  core-holder,  in  Figure  1 1 1 - 5 . 

The  microwaves  were  attenuated  in  accordance  with  Eq . ( 1 1 1 -8 ) 
and  were  received  on  the  opposite  side  by  the  horn  and  lens  antenna. 

A  thermistor  mount  (HP  R486A)was  used  for  detection  and  conversion  to 
a  signal  suitable  for  the  power  meter  (HP  431 C ) .  The  power  meter 
indicated  the  power  magnitude  in  milli-watts  and  this  power,  converted 
to  a  proportional  d.c.  voltage,  wasted  to  the  Y-axis  of  the  X-Y  recorder. 
A  close  up  of  the  microwave  instrumentation  is  shown  in  Figure  1 1 1-4 - 

111. 2. 3  Bench  Experiment 

A  semi-log  plot  of  attenuation  versus  power  transmitted  and 
received  by  the  horns  is  shown  in  Figure  1 1 1-7.  Power  received  is 
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FIGURE  I I 1-4 

CLOSER  VIEW  OF  MICROWAVE  INSTRUMENTATION 


FIGURE  III- 5 

HORN  AND  LENS  ANTENNAS  IN  POSITION 
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shown  along  two  straight  lines,  one  with  the  lenses  on,  and  one  without 
them.  This  bench  experiment,  with  the  antennas  three  inches  apart  in 
the  air,  clearly  demonstrates  the  advantage  of  using  a  lens  in 
conjunction  with  the  horn.  An  eighteen  percent  gain  in  received 
power  is  observed  because  the  spherical  wave  is  converted  to  a  plane 
wave,  thus  reducing  the  leakage  to  atmosphere. 

III. 3  FLUID  INJECTION  SYSTEM,  CORE  HOLDER  AND  ALLIED  EQUIPMENT 

A  Ruska  pump  was  used  to  inject  mercury  into  the  cylinders 
(Fig. 1 1 1-6 ) .  This  mercury,  in  turn,  displaced  either  oil  or  water, 
whichever  the  case  may  be.  This  liquid  was  injected  into  the  core 
holder  through  flexible  Tygon  tubing  and  the  usual  valves  and  fittings. 
Pressure  was  measured  just  before  the  inlet  to  the  core  holder  with  a 
Heise  pressure  gauge  having  a  precision  of  0.2  osi  (1.379  kPa). 

The  core  holder  had  a  rectangular  cross-section  so  that  the 
microwaves  would  traverse  the  same  thickness  along  the  height  of  the 
core.  The  core  holder  was  fabricated  from  PVC  and  no  appreciable 
attenuation  of  microwave  power  was  observed  through  it.  The  PVC 
sheet  was  1.27  cm  thick.  Two  core  holders  were  used  in  this  study. 
These  were  approximately  of  the  same  dimensions.  Details  concerning 
the  core  holders  are  available  in  Table  1 1 1 - 1 . 

TABLE  III-l 


Core  Holder 
No. 

Thickness 

(cm) 

Cross-sectional 
Area  (cm2) 

Length 

(cm) 

Run  Nos. 

1 

1.350 

6.859 

91.44 

8  to  1 5 

2 

1.270 

6.477 

91.93 

16,18-21, 

23,  24 
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The  liquids  were  collected  in  50  cc  tubes  which  were  placed 
in  a  fractional  collector.  The  core-holder  was  mounted  on  brackets 
with  wheels  and,  similarly,  the  fractional  collector  was  also  mobile. 
Both  these  units  were  driven  on  rails  by  a  chain  drive.  The  power  to 
the  chain  drive  was  imparted  by  a  variable  speed  d.c.  motor,  of 
one-sixth  horse  power,  through  a  series  of  sprockets.  A  limit  switch 
was  placed  along  the  rails  and  its  position  could  be  adjusted.  The 
limit  switch  was  used  to  shut  the  motor  off  when  it  was  pushed  by  one 
of  the  wheels.  A  photograph  of  the  complete  set-up  is  shown  in 
Figure  III-3. 

III. 4  A  TYPICAL  RUN 

III. 4.1  Preparation  of  the  Core 

Ottawa  sand  (80  -  120  mesh)  was  used  to  prepare  an 
unconsolidated  sand  pack.  The  mass  of  the  core  holder  was  determined 
before  and  after  it  was  packed.  The  core  holder  was  vibrated  using 
low  amplitude  vibrators  for  approximately  eight  hours.  During  this 
time  some  settlement  was  observed.  The  sand  was  then  topped-off  and 
the  inlet  end  plate  put  on. 

The  design  of  the  core  holder  was  such  that  the  outlet  end 
plate  was  an  integral  part  of  the  core  holder.  Sintered  teflon 
inserts  were  used  to  contain  the  unconsolidated  sand,  both  at  the 
inlet  and  outlet  ends.  The  thickness  of  these  inserts  was  0.3  cm  and 
they  were  highly  permeable.  The  sand  pack  and  the  inserts  were 
changed  for  every  run.  The  total  length  of  the  core-holder  was 
approximately  10  cm  longer  than  the  effective  length  (the  length 
containing  sand).  These  10  cm  were  at  the  inlet  end  and  contained 
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the  inlet  end  extension  (made  of  PVC)  plus  the  inlet  end  teflon  insert. 
This  measure  was  taken  to  keep  the  stainless  steel  bolts,  used  for 
securing  the  inlet  end  plate,  away  from  the  microwave  beam,  so  as  not 
to  cause  interference.  The  core  holder,  with  the  sand  in  it,  was  first 
weighed  and  then  mounted  on  the  brackets. 

The  position  of  the  limit  switch  and  the  speed  of  the  d.c. 
motor  were  adjusted  so  that  the  travel  time  of  the  core-holder  was 
15  seconds,  from  the  inlet  to  the  outlet  end,  relative  to  the  position 
of  horn  antennas.  This  was  done  because  the  speed  of  the  stylus  along 
the  X-axis  of  the  X-Y  recorder  was  kept  constant  at  1  inch/sec,  while 
the  length  of  the  X-axis  was  15  inches.  This  synchronized  the 
movement  of  the  core  holder,  relative  to  the  microwave  beam,  with  the 
travel  time  of  the  stylus  of  the  X-Y  recorder  in  the  x-direction, 
during  a  scan. 

The  microwave  equipment  was  turned  on,  and  the  horns  properly 

aligned  (Fig.III-5).  First,  a  scan,  with  only  air  between  the  horns, 

was  taken  and  the  power  meter  setting,  R  ,  and  X-Y  recorder  scaling 

factor,  S  ,  noted.  Then  a  dry  scan  was  carried  out  (S  =  0%)  along 
c  w 

the  length  of  the  core  holder  (Fig. I I 1-8) .  The  scan  was  conducted  by 
simultaneously  starting  the  motor  and  the  X-Y  recorder.  A  microwave 
response  was  obtained.  This  enabled  the  calculation  of  the  power, 

RDRY,  coming  out  of  the  core  at  0%  water  saturation. 

In  order  to  saturate  the  core  with  water  ora  glycerine-water 
mixture  (displacing  liquid  in  the  subsequent  waterflood),  a  vacuum  was 
pulled  on  the  core.  Water  was  injected  at  a  rate  of  1000  cc/hr  and 
the  volume  was  noted. 

A  preliminary  estimate  of  the  pore  volume  was  calculated 


using  the  equation, 
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Pore  Volume 


(III-l 9) 


where 


m 


mass  of  the  sand  in  the  core-holder,  gm, 


V  =  volume  of  the  core-holder,  cc. 


and 


p 


g 


grain  density  of  Ottawa  sand  =  2.64  gm/cc. 


When  this  calculated  pore  volume,  PV,  was  approached  during 


saturation  of  the  core,  the  pressure  rose  and  the  outlet  valve  was 
opened.  The  liquid  through-put  was  allowed  to  stabilize  and  the  total 
pressure  drop  across  the  core,  AP,  was  recorded  for  subsequent 
computation  of  the  absolute  permeability.  After  a  careful  material 
balance  on  the  liquid,  the  pore  volume,  PV,  was  determined. 


Because  of  equipment  considerations,  described  later  on, 


a  microwave  response  was  not  obtained  at  100%  S  when  water  was  the 

w 

saturating  liquid.  But  for  favourable  (y  /y  <  1)  water-floods  a 

mixture  of  glycerine  and  water  was  used  as  the  displacing  liquid. 

The  glycerine-water  mixture,  as  can  be  seen  later,  has  a  lower 

extinction  coefficient,  k,  and  consequently,  with  the  equipment  used, 

a  response  was  obtained  at  100%  S  .  (Figure  1 1 1 - 1 0 ) . 

w 

The  pack  was  kept  saturated  with  water  (or  glycerine-water 
mixture)  for  ten  to  fifteen  hours  so  as  to  render  the  sand  water-wet. 

m.4.2  The  Oil  flood 

In  order  to  establish  an  initial  water  saturation,  S  . , 

W  I 

an  oilflood  was  conducted.  The  rate  during  each  oilflood  was  the  same 


. 

. 
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in  a  particular  series  of  runs  (same  yQ/yw),  irrespective  of  the  rates 

used  in  the  subsequent  waterfloods.  This  was  considered  important  so 

as  to  have  similar  conditions  for  every  run  (in  a  series)  during  the 

establ ishment  of  S  . . 

wi 

Several  random  microwave  scans  were  conducted  during  an 
oilflood  (Fig.III-8).  This  helped  to  keep  a  check  on  the  displacement 
and  to  follow  the  flood  front  during  the  course  of  the  oilflood.  The 
oilflood  was  continued  till  no  water  was  produced.  At  this  time  a 
microwave  scan  gave  the  response  at  S  . . 

III. 4. 3  Power  Meter  Settings  and  Scale  Factors 

Before  describing  a  typical  waterflood,  it  is  necessary  to 

briefly  discuss  power  meter  settings,  X-Y  recorder  scale  factors  and 

the  various  techniques  used  in  getting  complete  microwave  responses. 

Figure  1 1 1 -7  indicates  that,  with  zero  attenuation  and  a 

lens-horn  antenna  combination,  the  power  measured  at  the  receiving 

horn,  through  air,  is  approximately  3.85  mW.  In  Appendix  C  an  analysis 

using  Eq . I I 1-7  is  undertaken  which  shows  that  the  attenuated  signal 

coming  out  from  the  core  can  vary  from  3  mW  at  0%  Sw  to  0.6(10  )  mW 

at  100%  S  ,  when  water  is  used  as  the  displacing  liquid, 
w 

During  a  single  scan  the  power  meter  setting,  R  ,  and  the 

X-Y  recorder  scale  factor,  S  ,  have  to  be  fixed.  The  microwave  power 

output,  which  ranges  from  3mW  to  0.6(10  )  mW  (say),  cannot  be 

recorded  with  a  single  combination  of  R  and  S  .  As  water  saturations 

vary  from  1  -  S  to  S.  during  a  waterflood,  such  a  spread  in 

o  r  wi 

microwave  responses  is  always  to  be  expected. 
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FIGURE  III-8  FACSIMILE  OF  RAW  DATA 
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Consider  a  power  of  R  mW  coming  out  of  the  core.  With  a 
power  meter  setting  of  R$  the  corresponding  d.c.  voltage  at  the  output 
of  the  power  meter  will  be,  (Fig.III-11), 

V(v0lts)  =  r^Svt  (nl-  20> 

where 

0  <  V  <  1.0 

This  d.c.  voltage  is  recorded  as  a  proportional  length  (in  inches) 

on  the  Y-axis  of  the  X-Y  recorder.  The  scale  factor,  S  ,  is  related, 

c 

to  the  voltage,  V,  and  the  distance  travelled,  Y,  by  the  stylus  of 
the  X-Y  recorder  along  the  ordinate,  through  the  equation. 

Y(inches)  =  sTmV/inchj  [Ref-  Fi9-HI-12]  (HI-21) 

where 

0  <  Y  <10  (because  the  length  of  Y-axis  is  10  inches) 
Re-writing  Eq . ( 1 1 1-21 ) , 

(S  ) ( Y ) 

v  =  ■  (II 1-22 ) 

10J 

Substituting  for  V  in  Eq . (III-20)  from  Eq . (III-22)  and  solving 
for  R  one  obtains 

R  =  (Rs)(Sc)(Y)10"3mW  (III-23) 

In  the  experiments  conducted  for  this  study,  power  meter 
settings,  R  ,  took  the  values  of  0.01,  0.1,  1.0  and  10  mW/V  while 
scale  factors,  S  ,  of  10,  50  and  100  mV/inch  were  employed. 

Considering  the  extrame  case  of  Sc  =  10  mV/inch  and  R$  =  0.01  mW/V, 
with  Y  being  read,  quite  accurately,  at  a  minimum  of  (say)  0.5  inches, 
then  the  value  of  R  measurable  using  this  equipment  can  be  computed 
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as  follows: 

Substituting  for  Sc,  Rg  and  Y  in  Eq. (III-23),  there  results, 

R  =  (0.01 )(10)(0.5)10"3mW 

or 

R  =  0.5(10'4)  mW 
or 

log10R  =  -4.3  (III- 24 ) 

A  cursory  glance  at  the  saturation  profiles  in  Appendix  A, 
for  those  runs  in  which  water  was  the  displacing  liquid,  will  show 
that  the  minimum  1  og-j  QR  response  is  close  to  the  number  computed  in 
Eq. (II 1-24 ) . 

III. 4. 4  The  Waterflood 

Before  a  constant  rate  waterflood  is  started,  the  total 
pressure  drop  at  the  initial  water  saturation  is  recorded.  This 
pressure  drop  is  required  to  compute  the  effective  permeabi 1 i ty ,  KQr, 
at  S  .. .  The  initial  water  saturation,  S  . ,  is  an  average  value  based 
on  a  knowledge  of  the  pore  volume,  PV,  and  the  material  balance  on 
the  oil  in  the  preceding  oil  flood. 

The  waterflood  was  now  started.  A  TIME-IT*  precision  watch 
was  used  to  record  the  time  from  the  inception  of  the  flood.  This 
gave  the  time  in  minutes  with  a  precision  of  one-hundredth  of  a 
minute. 

Usually  one  pore  volume  of  displacing  liquid  was  injected 
before  the  constant  rate  waterflood  was  culminated.  This  also  gave 


*  Precision  Scientific  Co.,  Chicago,  Ill.,  U.S.A. 
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the  opportunity  to  record  the  oil  production  after  the  breakthrough 
of  the  displacing  liquid. 

Saturation  scans  were  undertaken  at  suitable  time  intervals. 
To  obtain  a  complete  microwave  response,  at  any  instant,  the  core  was 
scanned  thrice.  The  three  scans  were  conducted  using  different 
power  meter  settings,  R  ,  but  usually  the  same  scale  factor,  S  . 

Figure  III-  9  is  a  facsimile  of  a  raw  data  sheet  as  observed  on  the 
X-Y  recorder,  at  forty  minutes  since  the  inception  of  the  waterflood, 
during  Run  18.  The  figure  shows  three  responses  at  power  meter 
settings  of  0.01,  0.1  and  1  mW/V,  the  scale  factor,  S  ,  being  the  same 
in  all  the  three  scans  (50  mV/inch).  The  importance  of  conducting 
three  scans  becomes  evident  upon  inspecting  the  raw  data.  For 
example,  the  microwave  response  at  an  R$  of  1  mW/V  is  virtually  zero 
for  the  region,  0.4  <  £  <  0.6,  but  in  the  same  region,  at  an  R$  of 
0.1  mW/V,  the  response  indicates  a  saturation  gradient  behind  the 
front.  This  saturation  gradient  could  not  have  been  recorded,  if 
different  settings  of  the  power  meter  were  not  used. 

The  total  pressure  drop  across  the  core  and  the  cumulative 
oil  or  oil/water  produced  were  also  recorded,  whenever  saturation 
scans  were  undertaken.  The  breakthrough  point  was  taken  as  the 
instant  when  the  first  water  droplet  was  seen  coming  out  of  the  core. 


III. 5  ANALYSIS  OF  ERRORS 

Errors  which  could  have  occurred  during  the  course  of 
experimentation  are  itemized  below.  Where  possible,  an  attempt  has 
been  made  to  estimate  a  quantitative  upper  bound  for  them. 
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FIGURE  III-9  FACSIMILE  OF  RAW  DATA 
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FIGURE  III-10  FACSIMILE  OF  RAW  DATA 
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Inches  on  Y-axis  R 
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POWER  METER 


</)  o  r-  o  (Settings  and  Input/Output) 


OUTPUT  (volts  d.c. ) 


FIGURE  III-ll 


X-Y  RECORDER 


(Scale  Factors  and  Input/Output) 


FIGURE  III-12 
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i)  During  a  water  flood,  for  a  particular  time,  three 
microwave  scans  were  conducted.  In  the  conversion  to  a  saturation 
profile  these  scans  were  assumed  to  depict  saturations  at  the  same 
instant,  which  was  not  exactly  true.  It  was  estimated  that  it  takes 
a  maximum  of  ninety  seconds  to  scan  three  times.  For  a  300  cc/hr 
waterflood,  using  core  holder  number  2  (Table  1 1 1 -1 ) ,  the  maximum 
percentage  relative  error,  with  respect  to  the  length  of  the  core, 
in  the  location  of  the  last  measured  saturation,  was  1.26  percent. 

ii)  Saturation  scans  appear  to  be  continuous  but,  in 
actual  fact  over  an  infinitesimal  period  of  time,  an  average  saturation 
in  an  elemental  volume  was  being  recorded.  The  horn  antenna  had  a 
dimension  of  3.462  cm  (Fig.  C-l )  in  the  vertical  plane.  The  vertical 
dimension  of  the  core  holder  was  5.08  cm.  Thus,  68%  of  the  height 
was  penetrated  and  was  taken  as  representati ve  of  the  elemental  volume 
computed  using  the  full  height  of  the  core.  Based  on  the  dimensions 
of  the  core  holder  number  2,  the  percentage  of  total  volume  which  was 
penetrated,  at  a  particular  instant,  was  computed  as  2.33%. 

Also,  for  a  core  holder  having  a  length  of  91.93  cm 
and  a  beam  width  of  2.15  cm  (Fig.  C-l),  forty  three  electrodes  would 
be  required  to  obtain  a  resolution  similar  to  microwave  scans,  if  the 
resistivity  technique  were  used. 

iii)  Energy  loss  through  attenuation  was  dissipated  in 
the  form  of  heat.  A  quantitative  analysis  carried  out  for  the  power 
levels  used  in  this  study  gave  a  temperature  rise  of  0.432°C/min. 

The  exposure  to  microwaves  for  the  whole  length  of  the  core  was 
fifteen  seconds  in  any  one  scan,  so  the  temperature  rise  was  not  an 
influence  of  any  consequence  [Parsons  ,  p . 305 J . 
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iv)  Another  source  of  error  could  have  been  the 
disturbance  of  the  horn  alignment  while  the  run  was  being  conducted. 
The  horns  were  aligned,  both  in  the  vertical  and  horizontal  planes, 
to  give  maximum  power  at  the  receiving  end  through  air.  This  had  to 
be  done  before  the  dry  core  scan  and  the  proper  alignment  should  have 
remained  till  the  last  water  flood  scan.  Utmost  care  was  taken  in 
this  regard,  and  periodic  checks  of  the  power  received  through  the 
air  were  carried  out.  Power  meter  readings  could  have  been  in  error, 
due  to  standing  waves,  if  the  load  (core)  was  not  matched.  A  tuner 
was  included  in  the  circuit  to  match  the  load. 

v)  The  displacement  of  oil  by  the  displacing  liquid 
usually  commenced  in  the  inlet  end  fitting  and  along  the  PVC  inlet 
insert.  These  have  a  cumulative  volume  of  7  cc.  This  volume  of 
displacing  liquid  was  injected  and,  as  a  consequence,  the  same 
volume  of  oil  produced,  before  the  actual  displacement  in  the  porous 
medium  commenced  at  £  =  0.  Due  to  this,  microwave  scans  measured 
saturation  profiles  which  did  not  measure  the  first  7  cc  of  volume 
injected,  although  this  volume  was  included  in  the  material  balance 
on  oil.  The  material  balance  was  corrected  when  experimental 
saturation  profiles  were  compared  with  the  numerical  results. 

Corrected  material  balance  values  were  used  in  the  data  bases 
compiled  for  the  numerical  simulators. 

vi)  At  the  moment  when  a  waterflood  commenced,  the  whole 
system  from  the  Ruska  pump  to  the  outlet  end  of  the  core  was  at 
atmospheric  pressure.  The  displacement  took  place  above  atmospheric 
pressure.  The  liquids  in  the  core,  theoretical ly ,  were  considered 
incompressible,  but  the  liquid  system  as  a  whole  including  the  mercury 
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and  the  liquids  in  the  cylinders  had  to  be  slightly  compressed.  The 
system  compressibil ity  shows  up  in  the  oil  produced  versus  pore  volume 
injected  curve.  Consequently,  the  linear  portion  of  the  oil  produced, 
t,  versus  pore  volume  injected,  t  • ,  curve  does  not  always  pass  through 
the  origin  [Figures  D-9  to  D-12]. 

vii)  The  Heise  pressure  gauge  used  to  measure  the  pressure 
at  the  inlet  end  (taken  as  aP,  as  outlet  end  was  at  atmospheric 
pressure)  had  a  full  scale  reading  of  100  psi.  The  precision  was 
0.2  psi.  Thus  errors  may  have  been  introduced  at  pressures  of  1  psi 
and  below. 

viii)  The  water  breakthrough  point  was  taken  at  the 
instant  when  the  first  drop  of  water  was  visible  at  the  outlet  end. 

The  outlet  end  fitting  had  a  volume  of  5  cc  so  the  breakthrough  from 
the  porous  medium  occurred  5  cc  before  that.  As  this  error  was 
consistent  in  all  the  runs,  the  comparison  between  runs  was  done 
without  regard  to  it.  For  absolute  comparison  with  numerical  results 
the  breakthrough  point  has  to  be  corrected  accordingly.  If  due  to 
the  capillary  discontinuity  at  the  outlet  end,  the  displacing  liquid 
broke  through  after  building  up  to  a  maximum,  then  this  delay  was  also 
considered  when  runs  were  being  compared  among  themselves. 
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CHAPTER  IV 


CONVERSION  AND  PREPARATION  OF  DATA  FOR  ANALYSIS 


IV. 1  INTRODUCTION 

Fifteen  displacement  runs  were  conducted  in  this  study. 

These  were  conducted  on  nearly  identical,  unconsolidated  sand  packs 
and  were  grouped  according  to  the  viscosity  ratios  (y  /y  ).  The  water 
used  as  the  displacing  liquid  at  adverse  viscosity  ratios  was  of  IN 
salinity.  At  favourable  viscosity  ratios  a  mixture  of  water  and 
glycerine  was  used.  This  mixture  served  two  purposes:  firstly  it 
resulted  in  a  high  viscosity  for  the  displacing  liquid  and  secondly 
the  extinction  coefficient,  k,  was  reduced.  This  enabled  a  response 
at  100%  saturation  of  the  displacing  liquid  to  be  obtained. 

Table  IV-1  gives  a  brief  description  of  the  runs  conducted. 
Comprehensive  data  sets  on  individual  runs  are  available  in  Appendix  B. 

The  interfacial  tension  between  the  displaced  and  displacing 
liquids  was  determined  using  a  Spinning  Drop  Tensiometer.  The 
kinematic  viscosity  was  measured  using  a  Cannon  Fenske  Viscometer. 

The  density  was  determined  using  a  specific  gravity  bottle  and  a 
precision  balance. 


IV. 2  CONVERSION  OF  MICROWAVE  RESPONSE  TO  A  SATURATION  PROFILE 


The  relationship  which  governs  microwave  attenuation  and 
saturation  is  given  by  Eq . (II 1-8) .  Rearranging  this  equation,  one 
obtai ns 


4irf  hkcjjS^ 

log10Io  "  2.303c 


(IV-1) 
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TABLE  IV-1 


SUMMARY  OF  RUNS  CONDUCTED 


Run  No. 

u  /  u 

Displacing  Liquid 

Displaced  Liquid 

Rate 

(cc/hr) 

8 

9.34 

IN  Saline  Water 

Dow-Corni ng 

200 

10 

9.34 

IN  Sal ine  Water 

Dow-Corning 

100 

11 

9.34 

IN  Saline  Water 

Dow-Corni ng 

300 

12 

73.33 

IN  Saline  Water 

Paraffin  Oil 

25 

13 

73.33 

IN  Saline  Water 

Paraffin  Oil 

50 

14 

73.33 

IN  Saline  Water 

Paraffin  Oil 

80 

15 

73.33 

IN  Sal ine  Water 

Paraffin  Oil 

100 

16 

12.25 

IN  Sal ine  Water 

Mineral  Oil* 

100 

18 

12.25 

IN  Saline  Water 

Mineral  Oil 

100 

19 

12.25 

IN  Saline  Water 

Mineral  Oil 

100 

20 

12.25 

IN  Saline  Water 

Mineral  Oil 

200 

21 

12.25 

IN  Saline  Water 

Mineral  Oil 

300 

9 

0.266 

Water  +  Glycerine 

Keros ine 

100 

23 

0.184 

Water  +  Glycerine 

Keros ine 

100 

24 

0.184 

Water  +  Glycerine 

Keros ine 

50 

*  Mineral  oil  was  a  blend  of  paraffin  oil  and  heptane. 
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A  constant  power  of  I  is  being  radiated  towards  the  core  holder  at 
a  fixed  frequency,  f.  It  is  assumed  that  the  porosity  thickness 
product  is  essentially  uniform  along  the  length  of  the  core.  The 
extinction  coefficient,  k,  is  constant  because  during  a  run  the  same 
dielectric,  either  water  or  a  glycerine-water  mixture, is  being 
injected.  Given  these  assumptions,  a  linear  relationship  exists 
between  lognnI  and  S. 

The  raw  data  facsimile  of  Figure  1 1 1-8  is  analysed  and 
converted  first.  An  average  value  for  R  (or  I  )  from  the  dry  scan 

X 

response  was  determined.  This  was  an  arithmetic  average  along  the 
length  of  the  core.  The  power  measured,  R,  was  computed  using 
Ea . (Ill -23 ) .  The  logarithm  of  the  average  value  of  R  gave  RDRY. 

In  all  the  runs  conducted,  RDRY  varied  between  +1.0  and  0.  An  axis 
was  laid  out  on  the  ordinate  between  +1.0  and  -5.0,  decreasing  in  the 
upward  direction.  *  The  abscissa  represented  the  normalized  distance,?. 

For  those  runs  in  which  a  glycerine-water  mixture  was  used 
as  the  displacing  liquid,  the  microwave  response  at  100%  Sw  was  also 
obtained  (Figure  III-10).  An  arithmetic  average  of  log  R  was  taken 
which  was  denoted  by  RSAT. 

The  initial  water  saturation  scan  was  converted  next.  By 
picking  suitable  points  along  the  response  curve, R  was  obtained 
using  Eq . ( 1 1 1-23 ) .  The  logarithms  of  the  measured  power  were  plotted 
against  the  normalized  distance,  ?.  The  initial  water  saturation  need 
not  be  uniform  along  the  length  of  the  core  and,  in  reality,  it  varied. 

The  variation  was  slight  when  the  oil  flood  viscosity  ratio  was  favourable, 
but  significant  (Run  9,  Run  23  and  Run  24)  when  it  was  adverse.  An 
arithmetic  average  of  log  R  was  computed  for  the  S^.  just  plotted. 
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The  average  value  of  log  R  was  termed  RSWI  and  was  taken  as  being 
representative  of  the  average  S  .  obtained  by  material  balance. 

The  waterflood  microwave  responses  (Figure  III-9)  were  then 
converted  to  saturation  profiles.  Usually  three  scans  were  conducted 
to  obtain  one  profile.  Points  were  picked  on  these  three  responses 
according  to  the  degree  of  resolution  at  different  normalized  distances. 
Only  one  response  is  read  to  obtain  R,  at  every  £.  Eq . (Ill -23 ) 
enables  the  conversion  of  these  responses  to  R.  A  plot  of  log  R 
versus  5  gives  the  saturation  profile  for  that  time. 

The  microwave  scan  at  residual  oil  saturation,  S„  ,  was  only 

or 

obtained  for  runs  conducted  at  favourable  viscosity  ratios.  Like  the 
initial  water  saturation,  SQr  also  need  not  be  uniform.  An  average 
log  R  for  the  SQr  response  was  determined  and  this  number,  RSOR,  was 
taken  as  being  representative  of  the  SQr  computed  from  material 
balance  considerations. 

For  all  the  runs,  RDRY  and  RSWI  were  computed,  but  RSOR  and 
RSAT  could  only  be  measured  for  the  favourable  viscosity  ratio 
waterfloods.  Figure  IV-1  shows  the  plot  of  log  R  versus  S  .  Three 
straight  lines  are  drawn  for  illustration.  The  line  for  saline  water 
represents  all  the  runs  conducted  at  adverse  viscosity  ratios.  The 
remaining  two  are  for  the  glycerine-water  mixtures  used  at  favourable 
viscosity  ratios.  Similar  data  are  found  in  Table  IV-2. 

Saturation  profiles,  for  all  the  runs  plotted  and 
calibrated,  are  available  in  Appendix  A.  The  number  written  beside 
each  profile  is  the  cumulative  pore  volume  produced,  x,  up  to  the 
time  the  scans  were  conducted.  These  are  uncorrected  values  of  x. 

In  the  profiles  shown  for  runs  9,  23,  24  and  15  two  numbers  are 
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TABLE  IV-2 


MICROWAVE  CALIBRATION  DATA 


Run  No. 

RDRY# 

RSWI+ 

8 

0.500 

10.82 

-0.450 

10 

0.470 

12.60 

-0.450 

11 

0.480 

9.00 

-0.450 

12 

0.520 

9.10 

-0.310 

13 

0.440 

7.90 

-0.280 

14 

0.340 

7.90 

-0.300 

15 

0.440 

8.47 

-0.350 

16 

0.273 

10.65 

-0.600 

18 

0.327 

8.70 

-0.500 

19 

0.317 

6.10 

-0.575 

20 

0.300 

11.10 

-0.550 

21 

0.460 

10.18 

-0.600 

Run  No.. 

RDRY  RSAT 

cJ-S 

•i — 

RSWI 

Sor 

RS0R 

9 

0.510  -3, 

.625 

17.90 

-0.250 

16.34 

-3.1 

23 

0.273  -3. 

.153 

23.80 

-0.550 

12.66 

-2.7 

24 

0.474  -3. 

,219 

24.35 

-0.400 

16.30 

-2.6 

*  OSwi}/n 

=  9.38%  | 

•  1 

+  (eRSWI )/n 

=  -0.451 

i  These 

values  were  used  to 

draw  the  straight 

#  (xRDRY)/n 

=  0.41  J 

1  line 

for  sal i ne 

water  displacements 

in  Fig.] 
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written  beside  the  saturation  profiles.  The  lower  value  is  the 
planimetered  area  under  the  profile.  By  definition  this  number  should 
be  exactly  equal  to  the  oil  produced  in  pore  volumes,  x.  But  due  to 
the  experimental  errors,  mentioned  before,  the  area  under  the  profiles 
and  the  pore  volume  of  oil  produced  are  not  exactly  equal.  Moreover 
t  has  also  to  be  corrected  for  the  volume  of  the  inlet  end  fittings. 


IV. 3  ERROR  ANALYSIS  ON  SATURATION  PROFILES 

Relative  errors  for  the  material  balance  were  computed  for 
runs  9,  23,  24  and  15  using  Eq . ( I V-2 ) . 

(t  -  Area)100 

Relative  Error  =  — -  ( I  V-2 ) 

Tcorr 


where 

Tcorr  =  {T  '  w}  (IV' 

The  relative  errors  are  plotted  against  the  corrected 
values  of  pore  volumes  of  oil  produced,  Tcorr>  and  are  depicted  in 
Figure  I V-2 . 


IV. 4  PORE  VOLUMES  INJECTED,  PORE  VOLUMES  OF  OIL  PRODUCED 
AND  EFFECTIVE  VISCOSITIES 

To  estimate  relative  permeabilities  from  displacement  data 
[Johnson29  et  al . ,  1959],  the  following  curves  have  to  be  generated 

on 

[Jones  and  Roszelle  ,  1978]. 

i)  Effective  viscosities,  I"1,  versus  pore  volumes 
injected,  x. 

ii)  Pore  volumes  of  oil  produced,  x,  versus  pore 
volumes  injected,  x  • 


RELATIVE  ERROR  (%) 


80 


7(L 


A 


ERROR  ANALYSIS  ON  MATERIAL  BALANCE 


(FIGURE  IV-2 ) 
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Data  sets  for  x  (N  /P V)  versus  x.  are  tabulated  in  Tables  D-l 

r  * 

to  D-l 5.  Plots  have  also  been  generated  and  are  placed  in  Appendix  D, 
as  Figures  D-9  to  D-l 2. 

The  effective  viscosities  are  defined30  by  Eq.(IV-4), 


where 


and 


<«?>&-) 

b 


ub  ' q  'VAP 


( IV-4) 


=  viscosity  of  the  liquid  used  to  determine  relative 
permeabil ity, 

aP  =  pressure  across  the  core  during  a  waterflood: 

it  varies  with  x,  and  is  tabulated  in  Appendix  D, 
Tables  D-l  to  D-l 5 , 

qb  =  rate  at  which  absolute  permeability  was  determined, 
q  =  waterflood  rate  for  a  particular  run, 

AP^  =  pressure  difference  across  the  core  when  one  liquid 
is  flowing  at  a  rate  q^. 

Define 


yb 


( I V  -  5 ) 


where 


C.  =  constant  for  the  i-th  run. 

The  values  of  this  constant  for  all  the  fifteen  runs  are  given  in 
Table  IV-3.  To  obtain  A"1,  multiply  C  by  AP  in  kPa. 

Plots  of  I-1  versus  x.  are  placed  in  Appendix  D  as 


Figures  D-5  to  D-8. 


« 


*  •  " 
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TABLE  IV-3 

VALUES  OF  CONSTANT  C  TO  DETERMINE 
EFFECTIVE  VISCOSITIES 


Run  No.  C(mPa-sec/kPa) 


8 

0.163 

10 

0.291 

11 

0.111 

12 

1.161 

13 

0.691 

14 

0.462 

15 

0.345 

16 

0.367 

18 

0.333 

19 

0.382 

20 

0.179 

21 

0.121 

9 

0.065 

23 

0.052 

24 

0.112 
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All  the  plots  found  in  Appendix  D,  whether  ^  versus  x .  or 
x  versus  t  .  ,  are  grouped  according  to  the  viscosity  ratio  used  during 
a  particular  series  of  waterfloods. 

IV. 5  PRESSURE  DIFFERENTIAL  AND  PORE  VOLUME  INJECTED 

In  order  to  compare  AP  versus  behaviour  vis  a  vis 
experimental  saturation  profiles,  plots  of  AP  versus  r-j  were  generated. 
They  were  also  grouped  according  to  their  viscosity  ratios.  The 
plots  are  shown  in  Figures  D-l  to  D-4  in  Appendix  D. 


CHAPTER  V 


NUMERICAL  SOLUTIONS  AND  COMPARISON  WITH 
EXPERIMENTAL  SATURATION  PROFILES 


V.l  INTRODUCTION 

The  one-dimensional  immiscible  displacement  equation  has 

1  c 

been  solved  numerically  by  many  researchers  such  as  Douglas  ,  et  al. 
[1958],  McEwen38  [1959],  Hovannessian28  [1961]  and  Flock19,  et  al . 
[1977]  in  the  Eulerian  domain  and  by  Payers  and  Sheldon1^  [1959]  and 
Bentsen  [1977]  in  the  Lagrangian  domain.  As  experimentally  measured 
saturation  profiles  are  now  available  (Appendix  A},  it  is  possible 
to  compare  the  experimental  profiles  with  solutions  of  the  one¬ 
dimensional,  immiscible  displacement  equation.  This  comparison  was 
undertaken  for  the  following  reasons: 

i)  To  gain  more  confidence  in  the  experimental  technique, 

ii)  To  determine  which  numerical  solution,  Lagrangian  or 
Eulerian,  compares  more  favourably  with  the 
experimental  saturation  profiles, 

iii)  To  gain  further  insight  into  the  physics  of  immiscible 
displacement  in  porous  media. 

V. 2  NUMERICAL  SOLUTIONS 

The  one-dimensional  immiscible  displacement  equation  as 
nondimensional ized  by  Bentsen  [1976]  may  be  written  and  solved  in 
the  following  two  formulations. 
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V.2.1  Lagranqian  [Bentsen^,  1977] 

This  formulation  is  in  the  fractional  flow  domain, 

{ f  ( S* ,  x*)},  and  the  following  three  equations  are  solved 
successively,  using  their  finite  difference  analogs,  to  arrive  at  a 
saturation  profile,  S*(£,  t*). 

Consider  the  frontal  advance  formula  [Buckley  and  Leverett^] 


(iL_) 

'  9t*' 


S* 


(v-i) 


Taking  derivative  of  both  sides  of  Eq.(V-l)  and  integrating 
between  the  limits  zero  and  t*,  the  first  equation  in  the  solution 
procedure  is  obtained, 


(V-2) 


This  gives  an  estimate  of  (3£/3S*)t*  at  every  node,  on  the 

saturation  axis,  at  the  end  of  a  particular  time  step.  Of  course, 

an  initial  estimate  of  32f/3S*2  is  required  to  start  the  numerical 

procedure.  The  second  of  the  equations  is  the  fractional  flow 

35  3 

equation  [Leverett  ]  as  adapted  by  Bentsen  in  Eq.(II-7),  but  written 
as 


G(S*) 


NcC(S*) 
(— ) 


(V- 3) 


The  values  of  (3S/3S*)^*  Eq. (V-2)  at  every  node  are  used 

in  Eq.(V-3)  to  obtain  a  current  value  for  f  (S*,  t*)  which,  in  turn, 

results  in  the  evaluation  of  an  improved  value  for  32fw/3S*2. 

Equations  (V-2)  and  (V-3)  are  iterated  upon  till  convergence  is 

obtained  on  3 2f, ./ 3 S*2 . 

w 
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In  this  formulation  S*  is  a  dependent  variable  so  no 
iterations  are  required  on  the  G(S*)  and  C(S*)  functions. 

Finally  to  obtain  a  relationship  between  z  and  S*,  at  a 
particular  time,  the  third  equation  to  be  solved  successively  is 


c(s*.  T*) 


=  s(s*,  0)  + 


dx* 


( V-4 ) 


which  is  obtained  by  integrating  Eq.(V-l)  between  zero  and  t*. 

£ 

The  boundary  and  initial  conditions  used  by  Bentsen  in 
the  solution  of  Equations  (V-2)  and  (V- 3)  were 


fw(0.  T*)  =  o 

fw(l  .  T*)  -  1 


T*  >  0 


(V-5) 
(V- 6) 


fw(S*,  0)  =  1  0  <  S*  <  1  ( V- 7 ) 

In  contrast  to  all  other  numerical  solutions  published 
till  that  time  this  particular  Lagrangian  formulation  allowed  the 
inlet  saturation  to  build-up  with  time  to  its  maximum  value.  The 
Lagrangian  solution  of  Fayers  and  Sheldon^  did  not  give  the  time 
elapsed  for  a  particular  saturation  profile.  This  shortcoming  was 
also  overcome  in  this  particular  [Bentsen  ]  Lagrangian  solution. 


V.2.2  Eulerian  (see  Appendix  E  for  details} 

The  Eulerian  formulation  in  the  saturation  domain  can  be 

obtained  by  differentiating  Eq.(II-7)  with  respect  to  £,  and  then 

substituting  for  9f  /9£  in  Eq .  ( 1 1 - 6 )  resulting  in  Eq.(V-8),  which  is 

w 


{NcC(S*) 


dG(S*)  3S* 
dS*  a? 


9  S*  _  _9_ 

9  T*  ZZ 


(V-8) 


* 
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19 

A  form  of  Eq.(V-8)  was  solved  by  Flock  ,  et  al.  using  a 
semi -imp] icit  integration  scheme.  They  assumed  an  instantaneous  jump 
of  the  inlet  end  saturation  to  its  maximum  value.  It  could  not  be 
ascertained  how  the  non-linearities,  C(S*)  and  G(S*),  in  Eq.(V-8), 

were  handled  in  their  solution  procedure. 

In  this  study,  Eq.(V-8)  was  solved  using  its  finite 
difference  analog  and  an  implicit  scheme  in  the  time  derivative.  An 
attempt  was  also  made  to  let  the  saturation  build  up  to  its  maximum 
value  at  the  inlet  end.  Appendix  E  gives  details  regarding  the 
treatment  of  the  G(S*)  and  C(S*)  functions. 

The  boundary  and  initial  conditions  used  to  solve  Eq.(V-8) 

were, 

S*U,  0+)  =  F(c)  0  <  %  <  1  (V-9) 

S*(l,  t*)  =  0  0  <  x*  <  T*bt  (V-10) 

The  initial  saturation  at  the  inlet  end,  SI,  was  determined 
by  numerical  trials  and  then  allowed  to  build-up  with  the  help  of 
Eq.(V-8).  The  particular  finite  difference  analog  of  this  equation 
is  given  as  Eq . (E-l 8)  in  Appendix  E.  After  S*  has  built  up  to  one  the 
third  boundary  condition  becomes 

S*(0 ,  x*)  =  1  ( V — 1 1 ) 

V. 3  DATA  BASES  FOR  NUMERICAL  SIMULATORS 

The  numerical  simulators  based  on  the  Lagrangian  and 
Eulerian  formulations  both  required  identical  data  bases,  in  so  far 
as  the  physical  and  geometrical  properties  in  the  runs  were  concerned. 
The  experimentally  generated  saturation  profiles  in  Runs  9,  23  and  24 
were  selected  for  comparison,  because  they  were  representative  of 


■ 


' 
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those  obtained  in  runs  made  at  favourable  viscosity  ratios.  For 
adverse  viscosity  ratios,  Runs  11  and  18  were  chosen.  This  choice 
was  based  on  an  inspection  of  the  saturation  profiles  generated  in 
all  the  runs. 

V . 3 . 1  Data  Obtained  by  Direct  Measurement 

The  area  and  length  of  the  core  was  taken  as  the  inner 
dimensions  of  the  core  holder.  The  viscosities  of  the  displaced  and 
displacing  liquids  were  measured  in  the  laboratory  at  room  temperature. 
At  the  end  of  every  oil  flood,  when  S.  was  established,  the  effective 

W  I 

permeability,  KQr,  to  oil  at  S  .  was  evaluated.  The  material  balance 

derived  S'.,  which  is  an  average  value,  was  used  as  input  to  the 

numerical  simulators.  For  Runs  9,  23  and  24  the  residual  oil 

saturation,  S’  ,  and  the  effective  permeability  to  water,  Kwr,  at  SQr 

was  also  obtained  experimentally.  This  was  facilitated  by  a  quick 

cessation  of  oil  production  after  breakthrough  of  the  displacing 

liquid  (glycerine  +  water).  For  adverse  viscosity  ratios  the 

waterflood  was  continued  until  about  one  pore  volume  of  water  had  been 

injected.  As  oil  production  was  still  continuing  at  that  time,  SQr 

and  K  could  not  be  determined  experimentally, 
wr  r  J 

V.3.2  Estimates  of  Residual  Oil  Saturation  and  Effective 
Permeabil ity  for  Adverse  Viscosity  Ratio 

In  Appendix  D,  tables  of  data  are  provided  for  each  run 
giving  values  for  dimensionless  pore  volumes  injected,  t. ,  cumulative 
volume  of  oil  produced,  N  ,  and  pressure  drop  across  the  core,  AP. 

r 

Extrapolation  of  the  1/t.j  versus  AP  curve  gave  an  estimate 

of  AP  at  infinite  pore  volumes  injected,  because  l/xi  tends  to  zero 


. 
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EXTRAPOLATION  OF  AP  -  1/T  RELATIONSHIP 
TO  INFINITE  PORE  VOLUMES  INJECTED 
(data  from  Run  1 1 ) 


FIGURE  V-l 


y  qL 
AaP 

(10~3) (0. 08333) (91 .44) 

(6.859) (13. 1 )103 

8.48(10'8)cm1 2 
8.48  ym2 


1  _ i _ 1 _ i _ I _ _  J 

2  4  6  8  10  1  2 
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as  t.  tends  to  infinity.  At  infinite  pore  volumes  injected,  the 

residual  oil  saturation,  SQr,  was  deemed  to  have  been  achieved.  So 

the  calculated  value  of  effective  permeability,  using  the  extrapolated 

value  of  AP,  was  taken  as  an  estimate  of  K  .  This  was  done  for  all 

wr 

the  runs  conducted  at  adverse  viscosity  ratios  and  the  resulting  K 

wr 

values  are  given  in  Tables  B-l  to  B-3  in  Appendix  B.  For  conducting 

numerical  runs  for  Runs  11  and  18,  similarly  generated  K ^  values  were 

used.  This  procedure  is  illustrated  for  Run  11  in  Figure  V-l . 

In  order  to  estimate  residual  oil  saturations,  S  ,  for 

or 

Runs  11  and  18,  relative  permeability  curves  were  generated  using  the 
Jones  and  Roszelle^  graphical  method.  Plots  of  X~ (see  Table  VI-3} 
and  t  versus  t.  were  drawn.  Tangents  at  points  on  the  after¬ 
breakthrough  portion  of  these  plots  were  then  drawn.  The  intersection 
of  these  tangents  with  the  ordinates  gave  values  for  the  effective 
viscosities,  and  saturations,  S1(r,  at  the  outlet  end  of  the  core, 
as  a  function  of  t.  .  These  values  enabled  the  calculation  of  the 
Keo(Sw)/K  anc*  ^ew(Sw)/K  functions  respectively.  Plots  of  these 
functions  versus  S are  shown  in  Figures  (V-2)  and  (V-3)  for  Runs  11 
and  18.  Next  a  line  parallel  to  the  saturation  axis  was  drawn,  so  as 
to  intersect  the  ordinate  at  Kwr/K.  The  K0W/ K  curve  was  extrapolated 
to  intersect  this  line  at  P.  A  perpendicular  was  dropped  from  point 
P  to  the  abscissa  intersecting  it  at  Q.  Point  Q  was  taken  as  an 
estimate  of  (1  -  S*or). 

V.4  ESTIMATION  OF  PARAMETERS  IN  RELATIVE  PERMEABILITY  MODELS 

3 

The  relative  permeability  models  introduced  by  Bentsen 


were 


, 


. 


. 
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RUN  NUMBER  11 


0 


.2 


.4 


.6 


1 
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Krw(S*)  =  a(S*)S*2  (V-12) 

and 

Kro(S*)  =  e(S*)(l  -  S*)NK  ;  NK  =  1  (V-13) 

where 

a  +  b  (1  -  S*) 

*  Waw  -  S*  U-U) 

w 

and 

a  +  b  S* 

b(S*)  =  --0--  +°$3  (V— 1 5) 

0 

The  functions  a(S*)  and  $(S*)  modify  the  shape  of  the  normalized 

relative  permeability  curves  to  account  for  the  influence  of 

wettability  and  pore  size  distribution. 

To  ascertain  representative  values  for  a  ,  b  ,  a  ,  b  and 

w  w  o  o 

NK,  numerical  experimentation  was  carried  out  on  the  data  base  for 
Run  11.  The  experimental ly  measured  saturation  profile  at  xcorr  of 
0.24  {t*  =  0.375}  was  chosen  to  match  the  numerical  profile  generated 
using  a  set  of  the  aforementioned  parameters  in  the  relative 
permeability  model. 

The  parameters  varied  were  a  ,  bQ  and  NK.  The  parameters 

a  and  b  were  assigned  the  value  unity  as  this  value  gave  a  good 
w  w 

match  of  the  model  { Eq . (V-12)}, with  the  relative  permeability  curve 

30 

for  water,  generated  using  the  Jones  and  Roszelle  method  {Fig.V-4}. 

For  a  particular  set  of  values  chosen  for  a  ,  bQ  and  NK  the  ^e0 (Sw)/K 

curve  and  the  saturation  profile  were  drawn  (Figures  V-4  and  V-5}. 

Two  criteria  were  assigned  for  the  selection  of  a  particular 

choice  of  a  ,  b  and  NK.  These  were: 
o  o 


' 


. 
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i)  The  saturation  profile  generated  using  the  selected 
values  of  aQ,  bQ  and  NK  should  match  the  experimentally 
measured  saturation  profile. 

ii)  Experimentally  measured  breakthrough  recoveries  should 

agree  closely  with  the  numerical  results  calculated 

using  the  chosen  values  of  a  ,  b  and  NK. 

o  o 

Both  the  Lagrangian  and  Eulerian  numerical  solution 
procedures  were  tested.  Table  V-l  shows  the  calculated  breakthrough 
recoveries  for  different  values  of  parameters  used  in  the  relative 
permeability  models.  Experimental  breakthrough  recoveries  are  also 
1 isted  in  Tabl e  V-l . 

Eulerian  solutions  {for  NK  =  2,  a  =1  and  b  =1}  gave  much 

o  o 

lower  breakthrough  recoveries  than  the  experimental  recoveries. 

Because  of  numerical  dispersion  the  foot  of  the  saturation  profile 
arrived  at  the  outlet  end  earlier  than  dictated  by  the  Lagrangian 
solution  of  the  immiscible  displacement  equation  (Figures  V-8  to  V-l 0) . 

The  saturation  profile  generated  using  the  Lagrangian 
procedure  {for  NK  =  1 ,  aQ  =  2,  bQ  =  -1 }  gave  the  closest  match 
(Fig.V-5)  to  the  experimental  profile  and  for  the  same  set  of 
parameters  the  breakthrough  recovery  agreed  well  with  the  experimental 
recovery.  These  values  for  the  parameters  {NK  =  1 ,  aQ  =  2,  bQ  =  -1 } 
were  used  to  calculate  numerical  profiles  which  were  in  turn 
compared  with  the  experimentally  measured  saturation  profiles 
{Figs. CV-11 )  to  (V-l 5) }.  The  choice  of  the  Lagrangian  procedure  was 
also  in  agreement  with  Fayers  and  Sheldon's^  conclusion  that 
Lagrangian  techniques  should  be  used  to  calculate  accurate  front 
profiles. 


'  *’#  *  -•**  I 
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BREAKTHROUGH  RECOVERIES  OF  SOME  RUNS  FOR  DIFFERENT  PARAMETERS 
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EFFECT  OF  DIFFERENT  SET  OF  PARAMETERS  ON  K  /K  CURVE 
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V. 5  ESTIMATION  OF  CAPILLARY  PRESSURE  NORMALIZING  PARAMETER 

V.5.1  Estimate  Based  on  Capillary  Pressure  Curve  Published  by  Hagoort 

An  imbibition  capillary  pressure  curve  was  reported  by 
Hagoort  °  [1974,  p.70.  Fig. 6]  for  crushed  Pyrex  glass  of  sieve 
fraction  between  450  and  715  ym.  The  functional  relationship  for  the 
curve  was 


J(s*)  =  (V-16) 

S*u! 

where 

Pr/i<7? 

j(S*)  =  ,  (V-17) 


35 

is  the  Leverett  J-function. 

Equating  Eq . (V-16)  and  Eq . (V-17)  one  obtains 


p  =  0J4x  s*-0.4  (V-18) 

c  A Tu 

Using  the  data  for  Run  11  from  Table  B-l  and  substituting 

for  interfacial  tension,  y,  absolute  permeability,  K,  and  porosity,  <f>, 

in  consistent  units  one  obtains 

P  =  0.00482244  S*"0,4  atm  (V-19) 

c 

As  mentioned  in  Chapter  II  Eq. (11-16)  ,  the  following 

capillary  pressure  model  is  considered  to  be  applicable: 

tt  =  -&nS* 

c 

or 

P  =  -a£nS*  (V-20) 

c 

Values  of  P  were  obtained  for  different  values  of  S*  using 
c 

the  functional  relationship  given  by  Eq . (V-19) .  A  plot  of  £nS*  versus 

P  was  drawn  using  these  values.  This  was  the  best  straight  line  for 
c 


. 


. 
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these  data  (Fig.(V-6),  straight  line  A}.  As  is  evident  from  Eq.(V-20), 
the  slope  of  this  straight  line  is  the  capillary  pressure  normalizing 
parameter,  a.  For  Pyrex  glass  beads  a  value  of  a  was  estimated  as 
0.003125  atm. 


V.5.2  Estimate  based  on  Capillary  Pressure  Curve  Published  by  Leverett 

oc 

Leverett  [p.161,  Fig. 4]  reported  a  capillary  pressure  - 
saturation  relationship  for  clean  unconsolidated  sands. 

Using  Eq . ( V- 1 7 )  and  the  data  base  for  Run  11,  as  before, 
the  relationship  between  the  Leverett  J-function  and  capillary 
pressure,  P  ,  is 

PC  =  0.03446  J(Sw)  atm  (V-21) 

Values  for  the  Leverett  J-function  were  read-off  the 
aforementioned  capillary  pressure  saturation  curve  and  Table  V- 2  was 
generated  using  Eq . (V-21 )  and  the  normalized  saturation  definition. 


Eq. (11-15) 

S* 


TABLE  V-2 

VALUES  FROM  LEVERETT  CAPILLARY  PRESSURE  CURVE 


K  VJJ  y  C 

(read-off  (Eq. V-21 } 

Leverett 's  work}  atm 

Sw 

S* 

(Eq. 11-15} 

£nS* 

0.2625 

0.0090 

0.80 

0.3000 

0.0103 

0.70 

0.953 

-0.048 

0.3250 

0.0112 

0.60 

0.797 

-0.227 

0.3875 

0.0133 

0.40 

0.484 

-0.726 

0.5000 

0.0172 

0.20 

0.172 

-1.760 

■ 


ESTIMATE  FOR  CAPILLARY  PRESSURE  NORMALIZING  PARAMETER  (a) 
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Using  the  tabulated  values  in  Table  V- 2,  £nS*  is  plotted 
versus  P  Fig.(V-6),  straight  line  B  and  the  slope  of  the  resulting 
straight  line  gave  a  value  of  0.004  atm  for  the  capillary  pressure 
normalizing  parameter,  a. 

V.5.3  Numerical  Experimentation  in  the  Estimation  of  Capillary 
Pressure  Normalizing  Parameter 

Having  obtained  an  approximate  estimate  of  a, -it  was 
possible  to  do  a  similar  analysis  as  was  done  with  the  parameters  of 
the  relative  permeability  models.  Once  again  the  same  experimental 
saturation  profile  for  Run  11  was  chosen  for  comparison.  Three  values 
of  a  were  used  and  saturation  profiles  were  generated,  using  the 
Lagrangian  procedure,  with  Run  11  as  the  data  base.  Comparison  . 
between  these  numerically  generated  saturation  profiles  and  the 
experimental  profile  is  shown  in  Fig.(V-7).  The  criterion  used  for 
choosing  a  was  that  the  slope  of  the  calculated  front  should  match 
that  of  the  experimental  profile.  For  Run  11  a  value  of  0.02  atm  for 
a  gave  the  closest  match  to  the  experimental  profile. 

For  other  runs  requiring  a  value  of  a  in  their  data  bases, 
similar  numerical  experiments  were  conducted.  These  values  of  a  are 
given  in  Table  V- 3. 

TABLE  V-3 

ESTIMATES  OF  CAPILLARY  PRESSURE 
NORMALIZING  PARAMETER 


Number 

0 

(atmospheres) 

11 

0.020 

18 

0.005 

9 

0.007 

23 

0.007 

24 

0.007 

’ 
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The  difference  in  the  values  of  a  estimated  using  capillary 

23  35 

pressure  curves  published  by  Hagoort  and  Leverett  with  that 
obtained  by  numerical  experimentation,  for  Run  11,  may  be  indicative 
of  a  different  wettability  characteristic  for  the  oil-sand  system 
used  in  this  study.  The  Leverett  J-function,  Eq.(V-17),  as  defined 
in  the  analysis  for  this  work,  does  not  include  the  effect  of 
wettabil ity. 

V . 6  COMPARISON  OF  EULERIAN  AND  LAG RANG  I AN  SOLUTIONS  OF  THE 
IMMISCIBLE  DISPLACEMENT  EQUATION 

Fayers  and  Sheldon^  [p . 1 53 ,  Fig. 17]  compared  saturation 
profiles  calculated  using  Lagrangian  and  Eulerian  formulations  of  the 
immiscible  displacement  equation.  Their  comparison  qualitatively 
indicated  the  relative  position  of  saturations  using  the  two 
approaches.  As  their  Lagrangian  scheme  could  not  give  the  time 
elapsed  for  a  particular  saturation  profile,  it  was  not  quantitatively 
evident  whether  the  Eulerian  and  Lagrangian  profiles  have  the  same 
relative  positions  at  equal  times.  As  a  procedure  based  on  a 
Lagrangian  formulation  [Bentsen"]  which  gave  the  elapsed  time,  for  a 
particular  saturation  profile,  was  available  it  was  decided  to  carry 
out  a  comparison  of  this  sort  with  the  solutions  based  on  the 
Eulerian  formulation. 

Data  sets  from  Runs  11,  18  and  24  were  used.  The  results 
from  the  Eulerian  procedure  were  obtained  first.  As  the  Eulerian 
procedure,  developed  and  programmed  for  this  study,  calculated 
saturation  profiles  at  multiples  of  Ax*  (time  step  size),  a  specific 
number  of  them  were  generated  until  breakthrough.  For  a  material 


k  '  I 

. 


COMPARISON  OF  EULERIAN  AND  LAGRAN6IAN  SOLUTIONS 


78 


■x-x-xxxxxx-x  Lagrangian 


79 


1.4  0.5  0.5  0.7  0.8  0.9 

NORMALIZED  DISTANCE 


FIGURE  V-10 


80 


NORMALIZED  DISTANCE 


81 


balance  check  the  area  under  these  profiles  was  also  calculated. 

The  maximum  material  balance  error  was  of  the  order  of  two  percent. 

The  profiles  were  then  compared  and  the  results  are  shown  in 
Figure  V-8  to  V-10. 

All  the  Eulerian  saturation  profiles  show  a  less  steep  front 
and  a  reversal  in  curvature  at  the  foot  of  the  front.  This  resulted 
in  early  breakthrough  of  the  displacing  liquid  and  pessimistic 
breakthrough  recoveries  of  oil  as  compared  to  experimental  recoveries 
and  those  calculated  using  the  Lagrangian  formulation.  The  results 
of  this  comparison  also  agree  with  what  was  inferred  by  Fayers  and 
Shel  don^. 

Based  on  these  comparisons  and  the  numerical  experimentation 
done  for  the  estimation  of  parameters  in  the  relative  permeability 
and  capillary  pressure  models,  it  can  be  inferred  that  solutions 
based  on  the  Lagrangian  formulation  should  be  used,  when  more  realistic 
frontal  shapes  and  breakthrough  recoveries  are  desired. 

V. 7  COMPARISON  OF  EXPERIMENTAL  SATURATION  PROFILES  AND 
LAGRANGIAN  SOLUTIONS 

Data  bases  required  to  support  the  numerical  procedure  were 

prepared  for  Runs  9,  23,  24,  11  and  18.  For  favourable  viscosity 

ratios  (Runs  9,  23  and  24)  all  the  data  sets  were  determined 

experimentally,  but  for  adverse  viscosity  ratios  (Runs  11  and  18), 

the  residual  oil  saturation,  SQr,  and  effective  permeability  to  water, 

K  ,  at  S  were  determined  by  using  the  extrapolation  methods 
wr  or 

mentioned  earlier  in  this  chapter.  The  dimensionless  pore  volumes 
produced,  t,  were  corrected  using  Eq . ( IV-3 ) .  These  corrected  values 


' 


' 
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are  written  beside  the  saturation  profiles  shown  in  Figures  V-II  to 
V-15.  In  the  experimental  saturation  profiles  given  in  Appendix  A, 
uncorrected  values  of  dimensionless  pore  volumes  produced,  t,  are 
written  beside  the  saturation  profiles. 

The  relative  permeability  and  capillary  pressure  models 

34  23  35 

were  taken  from  the  literature  5  and  preliminary  estimates  *  for 
their  parameters  arrived  at.  Representative  values  for  the  parameters 
were  obtained  by  numerical  experimentation. 

In  Runs  9,  23  and  24  the  oilflood  (while  establishing  S  .} 

was  carried  out  at  an  adverse  viscosity  ratio.  A  saturation  gradient 

in  the  established  S  .  is  evident  in  all  three  runs  (Figs.  V- 1 1  to 

V- 1 3 } .  In  the  data  sets  used  to  generate  numerical  saturation 

profiles  an  average  initial  water  saturation,  S’. ,  determined  from 

material  balance  was  used.  The  numerical  procedure  considered  this 

value  of  S  .  as  constant  throughout  the  length  of  the  core.  This  was 

not  the  case  in  the  actual  experiments.  Although  the  average  absolute 

errors,  in  the  material  balance,  (Fig.IV-2)  for  Runs  9,  23  and  24  were 

3.74%,  13.55%  and  3.29%  respectively,  the  experimental  saturation 

profiles  seem  to  run  faster  than  the  numerical  profiles  at  the  later 

stages  of  the  waterflood.  This  can  be  attributed  to  the  increasing 

S  *  towards  the  outlet  end  of  the  core.  Otherwise,  the  slopes  of  the 
wi 

profiles,  and  the  particular  relative  permeability  and  capillary 
pressure  models  chosen  seem  to  prevail  in  conjunction  with  the 
immiscible  displacement  equation,  in  so  far  as  the  overall  physics 
of  the  displacement  is  concerned. 
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For  the  adverse  viscosity  ratios  in  Runs  11  and  18  the 

initial  water  saturation,  S. ,  is  more  uniform  and  the  experimental 

profiles  tend  to  agree  much  more  closely  with  the  numerical  results 

{Figs.  V-14  to  V-l 5} .  Again,  for  those  particular  runs  selected  for 

numerical  comparison,  the  immiscible  displacement  equation  seems  to 

apply.  It  may  be  worthwhile  to  mention  here  that  as  the  immiscible 

displacement  equation  is  solved  using  normalized  saturations,  S*,  they 

have  to  be  de-normal i zed  to  saturations,  S  .  This  is  done  by  using 

w 

Eq. (11-15).  As  a  consequence  the  shape  of  the  profiles  is  a  strong 
function  of  the  residual  oil  saturation,  SQr.  An  experimental 
determination  of  SQr  would  have  given  more  confidence  in  the  analysis 
of  immiscible  displacement  theory  undertaken  here. 
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CHAPTER  VI 


DISCUSSION  OF  RESULTS 


VI. 1  INTRODUCTION 

Displacements  may  be  either  stable,  or  unstable.  If  the 
displacement  is  unstable,  viscous  fingering  is  completely  dominating 
the  process  and,  as  a  consequence,  displacement  behaviour  cannot  be 
predicted  by  solving  the  displacement  equation.  Even  though  the 
displacement  is  stable,  it  may  not  be  fully  stabilized.  Displacements 
are  said  to  be  stabilized  when  the  frontal  region  of  the  saturation 
profile  has  attained  (or  nearly  attained)  a  fixed  configuration. 

That  is,  the  displacement  is  stabilized,  if  the  saturation  profile 
is  close  to  that  predicted  by  the  conventional  Buck! ey-Leverett^ 
theory. 

i 

VI. 2  ATTAINMENT  OF  THE  STABILIZED  ZONE 

53 

For  adverse  viscosity  ratios,  Warren  and  Calhoun  [1955] 
have  reported  on  several  displacement  tests.  They  indicated  that, 
when  total  pressure  drop,  aP,  is  plotted  versus  time,  three  distinct 
regions,  prior  to  breakthrough,  could  be  identified  on  the  curve: 
a  primary  phase  where  the  pressure  is  building  up  due  to  compressibil ity 
effects,  a  non-linear  secondary  phase  where  the  pressure  declines  but 
the  frontal  region  has  not  yet  attained  a  stabilized  configuration 
and,  finally,  a  linear  decline  in  pressure  indicating  that  a  stabilized 
zone  has  formed.  Such  plots  for  adverse  viscosity  ratios  are  shown 
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; 
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in  Figure  D-l  to  D-3.  Breakthrough  points  are  also  indicated  on  these 
plots. 

Figure  D-l  shows  the  plots  for  the  displacement  of  Dow 

Corning  oil  (y  /u  =  9.34).  Run  10  (100  cc/hr)  indicates  that  a 

0  w 

stabilized  zone  has  not  yet  formed.  But  the  plots  for  Run  8 
(200  cc/hr)  and  Run  11  (300  cc/hr)  are  such  that  a  straight  line  can 
be  drawn  through  the  pre-breakthrough  data,  suggesting  that 
stabilization  of  the  frontal  region  has  been  attained.  Analysing  the 
saturation  profiles  for  this  series  of  runs  (Figures  A-l  to  A-3) 
leads  to  similar  conclusions.  The  frontal  regions  in  the  profiles 
shown  for  Runs  8  and  11  (Figs.  A-l  and  A-3)  seem  to  have  a  fixed 
configuration  while  the  profiles  in  Run  10  (Fig.  A-2)  were  still 
trying  to  stabilize,  when  breakthrough  occurred. 

In  Figure  D-2,  AP  versus  t.  plots  are  shown  for  Paraffin 
oil  displacements(uo/uw  =  73.33)  conducted  at  different  rates.  Run  12 
(25  cc/hr)  and  Run  15  (100  cc/hr)  indicate  that  a  stabilized  zone  may 
have  been  formed  prior  to  breakthrough,  while  Run  13  (50  cc/hr)  and 
Run  14  (80  cc/hr)  do  not  depict  stabilized  behaviour.  The  saturation 
profiles  (Figs.  A-4  to  A-7}  for  this  series  of  runs,  due  to  a  larger 
adverse  viscosity  ratio,  are  obscured  due  to  perturbations  caused  by 
microscopic  heterogeneties  in  the  porous  medium.  Hence  no  conclusion 
can  be  drawn  regarding  the  attainment  of  a  stabilized  zone  by  only 
looking  at  the  saturation  profiles. 

For  the  Mineral  oil  di spl acements,  Figure  D-3  shows  the  plots 
of  AP  versus  t..  In  this  series  of  runs,  all  the  plots  indicate  that 
the  displacement  is  stabilized  prior  to  breakthrough.  The  saturation 
profiles  for  Mineral  oil  displacements  (Figs.  A-8  to  A-l 2 }  show  that 
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for  Runs  16,  18  and  19,  all  of  them  conducted  at  100  cc/hr,  the 

frontal  region  has  formed  a  fixed  configuration  prior  to  breakthrough. 

The  saturation  profiles  for  Run  20  in  Fig.  A-ll  have  slopes  less  steep 

than,  say,  the  profiles  in  Fig.  A-9  for  Run  18.  The  only  difference 

in  these  two  runs  is  the  rate.  Run  20  (200  cc/hr)  should  have  steeper 

profiles  than  Run  18  (100  cc/hr),  as  dictated  by  immiscible  displacement 

53 

theory.  But  this  was  not  observed.  As  Warren  and  Calhoun's  analysis 
of  aP  versus  t.  plots  was  based  on  the  conventional  immiscible 
displacement  theory,  which  is  violated  by  the  profiles  of  Run  20,  any 
deductions  for  Run  20  and  Run  21  (which  was  conducted  at  an  even  higher 
rate  of  300  cc/hr)  regarding  the  attainment  of  a  stabilized  zone  will 
be  without  any  theoretical  basis. 


As  is  evident  from  the  foregoing  analysis,  AP  versus  t. 


behaviour  was  not  a  very  reliable  test  for  determining  whether 
stabilization  had  been  attained.  On  the  contrary,  it  suggested 
stabilized  displacements  even  for  those  runs  which  did  not  conform  to 
the  immiscible  displacement  theory.  Thus,  a  more  robust  method  of 
attack  was  required  to  analyse  the  displacement  tests  conducted  in 
this  study. 

VI. 3  ANALYSIS  WITH  THE  HELP  OF  DIMENSIONLESS  GROUPS 


Bentsen  5  has  shown  that  two  scaling  groups,  the  mobility 


ratio,  M  ,  and  the  capillary  number,  N  ,  are  important.  They  are 
■  c 


M 


r 


(VI-1) 


N 
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(VI-2) 
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These  scaling  groups  pertain  if  the  displacement  is  dictated  by 
immiscible  displacement  theory,  and  hence  is  stable,  but  not 
necessarily  stabilized. 

1  q  qq  a  r 

For  unstable  displacements  [Chouke  ,  Hagoori  ],  Peters^0 
has  proposed  a  stability  criterion  which,  for  rectangular  co-ordinates, 


is 


(M  -  1) 


Vy  L2L2 
x  y 


C*yK  (L2  +  L2) 
y  wrv  x  y' 


-  7T 2 


(VI-3) 


When  Ep  is  greater  than  tt2,  the  displacement  is  unstable 
and  immiscible  displacement  equations  do  not  apply.  The  capillary 
constant,  C*,  is  unique  for  the  fluids  and  porous  medium  used.  In 
this  study  Runs  8  to  11  used  the  same  oil  {Dow  Corning},  though  of 
different  viscosity,  and  the  same  type  of  sand  {Ottawa  sand}  as  used 
by  Peters^6  in  his  experiments.  For  this  system  he  estimated  a  value 
of  306.25  for  C*,  in  connate  water  bearing  displacements. 

C 

Bentsen  has  suggested  that  a  plot  of  breakthrough  recovery 

versus  reciprocal  capillary  number  should  approach  asymptotical ly  the 

1 0 

Buckl ey-Leverett  steady  state  breakthrough  recovery.  For  those 
relative  permeability  models  which  seem  to  apply  in  the  runs  conducted 
in  this  study,  the  steady  state  breakthrough  recovery  can  be  derived 

3 

in  a  manner  similar  to  that  of  Bentsen  . 

Consider  the  non-capillary  fractional  flow  function 


Fv/S*> 


M  K 
r  rw 


M  K  +  K 
r  rw  ro 


(VI-4) 


and  the  relative  permeability  models 
Krw(S*)  =  S*2 

Kr0(S*)  =  (1  -  S*) 


(VI-5) 


and 


(VI -6 ) 
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Substitute  for  the  relative  permeability  models  from  Eqs . (V 1-5)  and 
(VI-6)  in  Eq . (VI-4 )  to  obtain 


M  S*2 

F  (s*)  =  - - - 

M  S*2  +  (1  -  S*) 


(VI-7) 


or 


dpw(s*) 

dS* 


MrS*(2  -  S*) 


{MrS*2  +  (1  -  S*)}2 


(VI-8) 
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By  virtue  of  Welge's  tangent  construction,  the  flood  front 
saturation  can  be  related  to  the  fractional  flow  ^unction  as. 


’dFw(S*) 

dS* 


Sf 


Fw<S*f> 

Sf 


(VI-9) 


Substitute  for 


( 


dF  (S*) 


w 


dS* 


S* 


and  FW(S|)  from  Eqs. (VI-7)  and 


(VI-8)  in  (VI-9)  and  solve  for  S|  to  obtain 

t 

1 


s*  = 


AT 

r 


>  1  ;  and  S*  =  1  ,  M  <1 
r  t  r 


(VI-10) 


The  average  saturation  at  breakthrough,  S*^,  can  also  be 
related  to  the  flood  front  saturation  as 


S* 

^abt 


s*f 

W 


(VI-11) 


The  average  saturation  at  breakthrough  is  numerically 
equal  to  the  dimensionless  pore  volumes  injected  at  breakthrough, 
t£..  Substitute  for  S*  and  F(S$)  in  Eq. (VI-11)  to  obtain 

Dt  T  W  T 


T  C  j.  =  S* 


.* 

bt 


abt 


T* 

Tbt 


=  1 


jJ-  -  1}  Mr  >  1 


Mr  5  1 


(VI-12) 
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Breakthrough  recoveries,  x£t  for  Run  11  (Mr  =  8.24)  and 
Run  18  =  9.21),  using  Eq. (VI-12),  are  57.54%  and  55.5%,  respectively. 

These  are  shown  as  horizontal  lines  in  Figure  VI-1.  A  numerical 
experiment  was  undertaken,  using  the  data  bases  for  Runsll  and  18,  in 

which  was  plotted  versus  1/NC.  Maximum  breakthrough  recoveries  can 

51 

only  be  obtained  when  the  flood  is  stabilized  {Terwilliger  et  a  1.},  that  is, 
when  steady  state  (Bentsen  }  is  approached  asymptotically.  The  arrows 
in  Fig. VI-1  indicate  the  points  on  the  curve  at  which  the  two  runs 
were  actually  conducted.  It  is  obvious  that  a  stabilized  zone  has 
not  been  attained,  contrary  to  the  analysis  of  AP  versus  x.  curves  in 
Figures  D-l  and  D-3. 

The  capillary  constant,  C*,  for  the  Dow  Corning  -  Ottawa  sand 

system  has  been  quantitatively  determined  by  Peters^  and  the 

asymptotic  approach  to  the  steady  state  recovery  has  been  illustrated 

in  Fig. (VI-1)  at  a  capillary  number,  Nc,  of  0.01. 

Limits  on  N  for  stable  and  stabilized  displacements  can  be 
c 

evaluated  for  this  system.  Re-writing  Eq . ( VI -3 )  in  terms  of  N  , 
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Apply  the  stability  criterion  on  E  for  rectangular  co-ordinates, 
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RECIPROCAL  CAPILLARY  NUMBER  VERSUS  BREAKTHROUGH  RECOVERIES 
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Nc  1  0.01  (VI-18) 

if  a  stabilized  configuration  is  to  be  achieved  within  the  confines 
of  the  core.  Substitute  the  relevant  values  of  the  parameters  in 
Eq. (VI-17)  from  the  data  set  for  Run  11  and  obtain  the  limit  on  N  for 
stable  displacement 

Nc  >  0.0418  (VI-19) 

Comparing  inequalities  (VI-18)  and  (VI-19)  it  is  evident 

t 

that  for  the  geometry,  fluids  (Dow  Corning)  and  porous  medium  (Ottawa 
sand)  used  in  this  study,  a  stabilized  saturation  configuration  cannot  be 
obtained  within  the  length  of  the  core. 

Table  VI-1  gives  the  values  of  the  dimensionless  groups, 
the  experimental  breakthrough  recoveries  and  the  steady  state 
breakthrough  recoveries  using  Eq. (VI-12)  (but  converting  them  to 
percent  of  initial  oil  in  place}  for  all  the  runs  conducted. 

For  Runs  8  to  11 ,  all  have  Ep  much  less  than  the  maximum  of  t t2, 
indicating  that  these  three  runs  were  stable.  The  capillary 
number  in  each  of  the  runs  is  more  than  0.01  suggesting  that 
stabilization  has  not  been  attained.  This  conclusion,  which  also 


' 
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TABLE  VI-1 

VALUES  OF  DIMENSIONLESS  GROUPS  AND  RECOVERIES 
FOR  THE  RUNS  CONDUCTED 


Run  No. 

M 

r 

Nc 

ep 

Recovery 
%  IOIP 

Steady 

State 

Recovery  . 

%  IOIP  Vpw 

8 

6.03 

0.1453 

1.47 

51.3 

72.8 

10 

9.00 

0.4140 

0.82 

39.6 

63.6 

9.34 

11 

8.24 

0.1527 

2.15 

38.5 

63.2 

12* ** 

16.80 

0.4930 

1.27 

18.5 

47.1 

13* 

40.16 

0.6167 

2.13 

27.9 

31 .5 

73.33 

14* 

43.90 

0.4930 

3.19 

23.1 

30.3 

15* 

28.60 

0.2248 

4.27 

24.4 

37.0 

16 

8.62 

0.1209 

2.09 

51.8 

63.2 

18 

9 . 21 

0.1113 

2.30 

48.6 

60.2 

19 

12.50 

0.1815 

2.01 

41.9 

51.8 

12.25 

20 

11.89 

0.0806 

4.28 

44.5 

55.8 

21 

10.73 

0.0467 

6.34 

45.5 

57.6  ' 

9 

0.3278 

0.0234 

** 

79.0 

80.0  ) 

23 

0.2920 

0.0228 

** 

83.3 

83.4 

Favourabl e 

24 

0.2710 

0.0424 

** 

78.4 

78.5  * 

*  a  for  these  series  of  runs  was  estimated  as  0.017  atm 

**  these  runs  were  unconditionally  stable. 
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applies  to  Runs  8  and  11,  is  in  contradiction  with  the  analysis  of 
AP  versus  t.  curve  where  it  was  inferred  that  a  stabilized  zone  had 
been  formed. 

The  magnitude  of  Ep  for  the  other  two  series  of  runs 

(Paraffin  and  Mineral  oils)  was  determined  using  the  same  value  of  C* 

46 

as  estimated  for  Dow  Corning  oil  -  Ottawa  sand  system  by  Peters 
This  allowed  the  calculation  of  relative  values  for  Ep  for  each  run. 
The  absolute  magnitudes  of  Ep  may  not  be  correct  because  an 
experimentally  determined  value  of  C*  was  not  available  for  the 
Paraffin  and  Mineral  oil  systems. 

For  the  Paraffin  oil  displacements  the  magnitude  of  the 
capillary  numbers  suggest  that  stabilization  was  not  attained,  but  the 
estimated  values  of  Ep  indicate  that  the  runs  were  stable.  A  similar 
analysis  suggests  that  the  same  conclusions  can  be  drawn  with  respect 
to  the  Mineral  oil  displacements  in  Runs  16,  18  and  19.  That  is,  they 
are  stable  but  stabilization  has  not  been  attained.  Runs  20  and  21 
have  a  lower  Nc  but  a  relatively  higher  value  for  Ep.  The  recoveries 
are  also  less  than  those  obtained  in  Runs  16  and  18.  This  suggests 
an  unstable  displacement.  The  saturation  profiles  (Figs.  A-ll  and 
A- 12}  further  substantiate  this  belief. 

Thus  it  is  possible  to  put  bounds  on  C*  for  the  Mineral  oil 
-  Ottawa  sand  system.  As  Run  18  (100  cc/hr)  is  stable  (see  saturation 
profiles  in  Fig.A-9),  Ep  should  be  less  than  tt2,  while  Run  20  is 
unstable  (Fig. A-ll)  indicating  E  should  be  greater  than  n2  which 

r 

suggests  that  the  capillary  constant,  C*,  must  lie  within  the  range 
given  below 


77.22  <  C*  <  141.45 


(VI-20) 
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VI. 4  EXPERIMENTAL  SATURATION  PROFILES  AND  BREAKTHROUGH  RECOVERIES 

The  stable  immiscible  displacement  theory  suggests  that  there 

should  be  an  increase  in  recovery  with  decreasing  capillary  number, 

for  the  same  mobility  ratio,  till  a  stabilized  displacement  is  attained. 

The  definition  of  the  capillary  number,  Eq . ( V I -2 )  shows  that  an 

increase  in  length  of  the  core  will  decrease  the  capillary  number,  N  . 

32 

Jones-Parra  et  al.  have  reported  that  when  they  increased  the 
length  of  the  core,  for  the  same  rate,  the  recovery  was  increased. 

The  capillary  number  is  also  a  function  of  the  capillary 
pressure  normalizing  parameter,  a,  and  effective  permeability  to  water 
at  residual  oil  saturation,  K  .  These  two  variables  were  not 

W  r 

measured  experimentally.  Their  estimates  were  obtained  by  the  procedure 
described  in  Chapter  V.  Similarly,  mobility  ratio  is  a  function  of 
^wr*  The  estimates  of  effective  permeability  at  residual  oil 
saturation,  Kwr,  suggested  that  it  was  influenced  by  rate.  As  a 
consequence  the  mobility  ratio  would  also  be  influenced  by  rate. 

It  may  also  be  that  the  Kwr  was  estimated  at  different  residual  oil 
saturations.  If  this  is  the  case  then  SQr  could  be  a  function  of  rate. 
Because  of  these  uncertainties  in  estimating  ,  its  value  should, 
ideally,  have  been  determined  experimentally. 

In  the  Dow  Corning  oil  displacements,  Run  8  (200  cc/hr)  gave 
higher  recoveries  as  compared  to  Run  10  (100  cc/hr)  and  Run  11 
(300  cc/hr).  A  lower  estimate  of  mobility  ratio  for  Run  8  explains 
the  higher  recovery  and  also  suggests  a  higher  floodfront  saturation, 

S|  { Eq. (VI-10)} .  This  difference  in  is  evident  when  the 
saturation  profiles  for  Run  8  (Fig.A-1)  and  Run  11  (Fig.A-3)  are 
compared.  An  outlet  end-effect  was  observed  in  Run  10.  The  recovery 
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would  have  been  much  lower  without  the  end-effect.  Run  10  has  a 
higher  Nc  and  hence  the  frontal  region  is  far  from  attaining  a 
stabilized  configeration.  This  is  in  contrast  to  Run  8  and  11  which 
are  closer  to  stabilization.  This  is  substantiated  by  the  saturation 
profiles  shown  in  Figs.  A-l  to  A-3. 

On  the  basis  of  Ep  values  {Table  VI-1}  the  Mineral  oil 
displacements  {Figs.  A-8  to  A-10),  Run  16,  Run  18  and  Run  19  are 
considered  to  be  stable.  Their  breakthrough  recoveries  are  consistent 
with  the  magnitudes  of  N  and  M  .  Runs  20  and  21  were  found  to  be 
unstable  with  a  high  Ep  and  lower  breakthrough  recoveries.  The 
significant  difference  in  the  slopes  of  the  saturation  profiles,  when 
Run  18  and  Run  20  are  compared  {Figs.  A-9  and  A-l 1 } ,  indicates  that 

a  different  mechanism  is  operating  in  each  run. 

6 

Bentsen's  contention  that  the  larger  the  value  of  Mr  the 
earlier  steady  state  is  approached  is  substantiated  when  comparing 
the  runs  for  the  Paraffin  oil  displacements.  Run  13  (50  cc/hr)  and 
Run  15  (100  cc/hr)  have  different  mobility  ratios  and  capillary 
numbers.  Run  15,  with  a  lower  N  should  be  closer  to  steady  state, 
but  is  not  (compare  the  experimental  and  steady  state  breakthrough 
recoveries  for  each  run),  because  of  the  lower  M  .  Run  12  and  Run  14 
have  exactly  the  same  N  but  Run  12  is  farther  away  from  steady  state 
than  Run  14,  because  Run  14  has  a  much  higher  M  . 

VI. 5  COMMENTS  ON  EXTERNAL  DRIVE  TECHNIQUES  FOR  THE  DETERMINATION 
OF  RELATIVE  PERMEABILITIES 

30 

The  graphical  techniques  proposed  by  Jones  and  Roszelle 
were  used  to  generate  the  relative  permeability-saturation  relation¬ 
ships  depicted  in  Chapter  V.  Relative  permeability  curves  were 


. 

. 
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determined  for  Runs  11  and  18.  The  basis  for  external  drive 
techniques  is  the  non-capillary  Buckley-Leverett^  theory.  To 
strictly  follow  the  theory  and  use  displacement  data  to  generate  KeQ/K 
and  K  ,/K  curves  the  displacement  should  have  achieved  steady  state 
prior  to  breakthrough.  That  is  to  say  that  N  should  have  been  less 
than  0.01,  which  was  not  the  case.  This  may  explain  the  scatter  in 
the  data  of  K  /K  versus  Sif  curve  {Figs.  V- 2  and  V-3 } . 

Nevertheless  the  graphical  techniques  were  used.  These 
techniques  require  that  effective  viscosity,  X~\  and  pore  volumes 
produced,  t,  be  plotted  versus  pore  volumes  injected,  t..  Moreover 
the  after  breakthrough  portions  of  these  curves  are  used  for  the 
determination  of  relative  permeabilities.  A  close  inspection  of  the 
-  t.  and  t  -  x.  curves  indicate  (Figs.  D-5  to  D-7  and  Figs  D-9  to 
D-ll}  that  there  is  no  substantial  difference  in  the  behaviour  of  these 
curves,  for  the  same  viscosity  ratio,  after  breakthrough  of  the 
displacing  liquid1.  The  only  exceptions  to  this  argument  are  the  Dow 
Corning  oil  displacements  {Figs.  D-5  and  D-9}  in  which  Runs  10  and  11 
are  similar  but  Run  8  is  different.  This  implies  that  essentially 
the  same  relative  permeability  curves  will  be  generated  for  a 
particular  viscosity  ratio.  But  when  the  saturation  profiles  of  any 
two  runs,  in  the  same  series,  say.  Run  18  and  Run  20  are  compared 
{Figs.  A-9  and  A-l 1 } ,  it  is  apparent  that  substantially  different 
relative  permeability  curves  would  be  needed  to  match  the  saturation 
profi 1 es . 

The  drawbacks  of  the  external  drive  technique  were  also 

1  Aq 

reported  by  Archer  [1973]  and  Sigmund  and  McCaffery '  [1977]. 
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Archer^  suggested  the  use  of  a  numerical  simulator  to  match  the  oil 

production  and  pressure  histories  with  a  suitable  choice  of  relative 

permeability  models.  His  contention  was  that  smooth  relative 

permeability  curves  will  be  obtained,  in  contrast  to  the  anomalous 

29 

ones  generated  using  the  JBN  method.  As  found  in  this  study,  even 
if  the  recoveries  and  pressure  histories  are  matched  for  a  particular 
set  of  relative  permeability  curves  this  does  not  imply  that  the 
saturation  profiles  will  also  be  matched  using  the  same  set  of 
relative  permeability  curves. 

VI. 6  ANALYSIS  OF  RUNS  AT  FAVOURABLE  MOBILITY  RATIOS 

The  capillary  numbers  for  Run  9,  Run  23  and  Run  24 
{Table  VI-1}  have  values  closer  to  0.01  implying  that  these 
displacement  runs  were  nearly  stabilized. 

Figure  D-4  shows  the  plots  of  pressure  versus  pore  volumes 
injected  for  each  of  the  three  runs  conducted  at  favourable  mobility 
ratios.  The  plots  do  exhibit  a  non-linear  secondary  phase  and  then 
a  linear  stabilized  zone,  prior  to  breakthrough.  The  saturation 
profiles  also  suggest  {Figs.  A- 1 3  to  A- 1 5 }  that  after  some  time  they 
do  attain  a  stabilized  configuration. 

Runs  23  and  24  have  essentially  the  same  mobility  ratio  but 
Run  24  has  twice  the  capillary  number  of  Run  23.  Both  are  nearly 
stabilized  but  the  slope  of  the  stabilized  front  at  the  lo wer  N  is 
steeper.  This  is  consistent  with  the  immiscible  displacement  theory 
which  suggests  that  as  capillary  forces  increase  {higher  Nc>  the 
frontal  region  becomes  more  diffuse. 


. 
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The  non-linear  secondary  phase  on  the  aP  versus  x.  plot, 
and  the  saturation  profiles  corresponding  to  that  phase,  are  in 
accordance  with  the  transient  behaviour  of  the  fractional  flow 
function  shown  in  Fig. I I - 4 . 

The  breakthrough  recoveries  in  Run  9,  Run  23  and  Run  24  are 
consistent  with  the  magnitudes  of  M  and  N  .  It  could  be  said  that 
these  runs,  conducted  at  favourable  mobility  ratio,  validate  the 

53 

arguments  for  aP  versus  x.  curve  put  forward  by  Warren  and  Calboun  , 

c 

and  substantiate  Bentsen's  contention  that  steady-state  solution  is 
an  acceptable  approximation  to  the  displacement  equation  at  low  N  . 

It  is  appropriate  to  mention  that  external  drive  techniques 
cannot  be  applied  here  as  oil  production  ceased  immediately  after 
breakthrough.  This  is  evident  from  T""*  and  x  versus  x-  curves  in 
Figs.  D-8  and  D-12,  respecti vely. 


CHAPTER  VII 


CONCLUSIONS  AND  RECOMMENDATIONS 

VII. 1  SUMMARY  AND  CONCLUSIONS 

An  attempt  has  been  undertaken  to  obtain  a  better 
understanding  of  the  immiscible  displacement  process  by  measuring 
saturations  as  a  function  of  time  and  distance  using  a  relatively  new 
technique  based  on  the  principle  of  microwave  attenuation.  The 
measured  saturation  profiles  were  compared  with  numerical  solutions 
of  a  Lagrangian  formulation  of  the  displacement  equation  to  yield  a 
quantitative  measure  of  the  theory. 

An  overall  analysis  in  the  light  of  published  scaling 
groups  and  a  stability  criterion  was  also  undertaken  on  saturation 
profiles,  pressure  plots  and  breakthrough  recoveries. 

The  experimental  and  numerical  results  of  this  study, 
together  with  the  analysis  undertaken,  led  to  the  following 
concl usions : 

1.  The  solution  to  the  Lagrangian  formulation  of  the  immiscible 
displacement  equation  compared  more  favourably  with  the 
experimental  saturation  profiles  than  did  the  solution  based  on 
the  Eulerian  formulation. 

2.  The  one-dimensional,  immiscible  displacement  theory 
unconditionally  represents  displacements  carried  out  at  a 
favourable  mobility  ratio.  For  adverse  mobility  ratio 
displacements  in  horizontal  systems,  a  stability  criterion  must 
be  satisfied  before  the  theory  can  be  said  to  represent  the 
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displacement  process. 

3.  A  stabilized  configuration  of  the  frontal  region  was  not  attained 
in  the  adverse  mobility  ratio  displacements  carried  out  in  this 
study. 

4.  In  the  Dow  Corning  oil  displacements,  a  stable  displacement 
together  with  a  stabilized  configuration  of  the  frontal  region 
was  not  possible,  for  the  core  dimensions  and  the  sand  used  in 
this  study.  That  is,  the  bounds  on  the  capillary  number,  N  , 

V»» 

for  stabilized  and  stable  displacements  were  mutually  exclusive. 

5.  The  capillary  constant,  C*,  was  found  to  lie  in  the  range 

77.22  <  C*  <  141.45 

for  the  Mineral  oil  -  Ottawa  sand  system. 

6.  The  capillary  pressure  normalizing  parameter  was  estimated  by  a 
history  match  of  the  saturation  profiles.  The  slope  of  the  front 
depended  on  the  magnitude  of  the  capillary  pressure  normalizing 
parameter,  a. 

7.  This  study  experimental ly  verified  that,  for  the  same  capillary 

number,  N  ,  the  difference  between  steady  state  and  transient 
c 

breakthrough  recoveries  is  higher  for  lower  mobility  ratios. 

8.  The  experimental  results  obtained  in  this  study  suggest  that  the 
displacement  process,  at  favourable  mobility  ratio,  can  be 
approximated  by  the  steady  state  solution  of  the  displacement 
equation,  provided  the  capillary  number  is  sufficiently  small. 

9.  The  scatter  in  the  data  for  the  oil  relative  permeability  curve, 
calculated  using  the  external  drive  technique,  may  be  due  to  the 
fact  that  the  runs  were  not  stabilized.  This  suggests  that  a 
stabilized  displacement  is  necessary,  if  the  external  drive 
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technique  for  the  determination  of  relative  permeabilities  is 
used. 

10.  It  may  be  difficult,  if  not  impossible  to  infer  the  correct 
relative  permeability  curves  from  effective  viscosity  and 
cumulative  oil  production  histories  by  applying  the  external  dri 
technique  used  in  this  study.  This  is  so  because,  while  the 
saturation  profiles  differed  substantially  in  a  given  series  of 
runs  (suggesting  different  relative  permeability  curves),  the 
effective  viscosity  and  cumulative  oil  production  histories  only 
differed  slightly  for  a  given  viscosity  ratio. 

VI I. 2  RECOMMENDATIONS 

The  following  suggestions  may  help  in  improving  upon  the 
results  of  this  study: 

1.  The  oscil lator/RF  unit  used  in  this  study  did  not  provide 
sufficient  power  to  measure  water  saturations  over  the  entire 
range.  It  is  recommended  that  a  higher  microwave  power  output 
unit,  such  as  a  Klystron,  be  used  in  future  studies. 

2.  Raw  microwave  data  were  converted  to  saturation  profiles  by 
manually  taking  the  logarithm  of  the  power  measured.  In  order 
to  reduce  the  manual  work  required  in  this  conversion,  it  is 
recommended  that  a  log  amplifier  be  placed  in  the  circuit. 

3.  An  extrapolation  technique  was  used  to  obtain  SQr  and  As 

there  is  some  question  as  to  the  accuracy  of  values  obtained  in 
this  way,  it  is  recoimended  that  SQr  and  Kwr  be  experimental ly 
determined  in  future  studies. 

4.  It  was  found  that  the  limits  on  the  capillary  number  for  a 
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stabilized  and  stable  displacement  were  mutually  exclusive  in 
some  of  the  runs  conducted  at  adverse  viscosity  ratios.  It  is 
recommended,  therefore,  that  some  care  be  taken  to  design  a  flood 
where  both  criteria  can  be  met,  if  external  drive  techniques  are 
used  to  estimate  relative  permeability  curves. 

5.  Runs  16,  18  and  19  were  replicates.  Run  19  was  considerably 
different  than  the  other  two.  While  this  difference  may  be 
attributed  to  the  difference  in  absolute  permeabi 1 i ty ,  it  may  also 
be  attributed  to  heterogeneous  wettability.  It  is  recommended 
that  the  sand  packs  used  in  future  studies  be  treated  so  as  to 
assure  homogeneous  wettability. 

6.  Parameters  used  in  the  capillary  pressure  and  relative 
permeability  models  were  obtained  by  matching  numerically  derived 
profiles  with  those  obtained  experimentally .  In  future  studies, 

it  is  recommended  that  these  parameters  be  estimated  independently. 

7.  The  microwave  technique  used  in  this  study  is  based  on  the  decrease 
in  amplitude  which  occurs  when  a  microwave  signal  passes  through 
the  core.  Saturations  can  also  be  estimated  by  measuring  phase 
shifts  in  the  microwave  signal.  It  is  recommended  that  an 
investigation  be  undertaken  with  a  view  to  ascertaining  which  of 
these  two  methods  is  better,  and  whether  they  may  be  used  together 
to  study  three-phase  flow,  for  example. 

8.  It  was  estimated  that  only  68%  of  the  sample  height  was  intercepted 
by  microwaves  in  the  vertical  plane,  which  may  not  be  true.  It 

is  recommended  that  an  experiment  be  carried  out  to  study  the 
variation  in  microwave  response  when  the  core  is  scanned  in  the 
vertical  plane.  This  experiment  may  suggest  that  improvements 
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in  the  design  of  the  horn-lens  antenna  are  necessary  if  a 
completely  representative  cross-section  of  the  core  is  to  be 
measured.  For  example,  a  larger  aperture  or  a  narrower  microwave 
beam  may  be  necessary. 

9.  Gravitational  segregation  may  be  a  problem  at  lower  displacement 
rates.  If  such  is  the  case,  it  is  recommended  that  the  smaller 
dimension  of  the  core  cross-section  be  placed  in  the  vertical 
plane.  As  a  consequence,  the  transmitting  and  receiving  horn 
antennas  have  to  be  placed  above  and  below  the  mobile  core-holder. 
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FIGURE  A-l  EXPERIMENTAL  SATURATION  PROFILES 
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FIGURE  A-2  EXPERIMENTAL  SATURATION  PROFILES 


Run  Number  1 1 
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FIGURE  A-3  EXPERIMENTAL  SATURATION  PROFILES 


Run  Number  12 


121 


MS  ‘uojieanies 

*• 

If)  CO  W  1- 

o  d  o  o  6  o0 


o  o  o  o  o 

■  ■  a  • 

^  CO  CM  r- 

i  •  •  i 


a  Boi 


FIGURE  A-4  EXPERIMENTAL  SATURATION  PROFILES 


Run  Number  13  “1 0.5 
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FIGURE  A-5  EXPERIMENTAL  SATURATION  PROFILES 


Run  Number  14  -(0.5 
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FIGURE  A-6  EXPERIMENTAL  SATURATION  PROFILES 


Run  Number  15 
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FIGURE  A-7  EXPERIMENTAL  SATURATION  PROFILES 
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FIGURE  A-8  EXPERIMENTAL  SATURATION  PROFILES 


Run  Number  18 
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FIGURE  A-9  EXPERIMENTAL  SATURATION  PROFILES 


Run  Number  19  nl.O 
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FIGURE  A- 1 0  EXPERIMENTAL  SATURATION  PROFILES 
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FIGURE  A-ll  EXPERIMENTAL  SATURATION  PROFILES 


Run  Number  21  -|1.0 
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FIGURE  A-l 2  EXPERIMENTAL  SATURATION  PROFILES 
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FIGURE  A— 1 3  EXPERIMENTAL  SATURATION  PROFILES 
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FIGURE  A- 1 4  EXPERIMENTAL  SATURATION  PROFILES 
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FIGURE  A-15  EXPERIMENTAL  SATURATION  PROFILES 
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C.l  DESIGN  OF  THE  LENS  ANTENNA 
C.  1 . 1  Basis 

The  design  of  the  lens  antenna  is  based  on  the  principle  of 
refraction  in  optics.  A  spherical  or  cylindrical  wave  has  to  be 
incident  on  the  excited  side  of  the  lens  and  emerges  into  free  space 
according  to  the  equation 


=  c_  =  3  CIO10) 

%  '  V 


n  =  =  refractive  index  f  1 

c  =  velocity  of  light  cm/sec 

v^  =  phase  velocity 

$  =  angle  of  incidence 

i p  =  angle  of  refraction 

In  this  work  the  radiator  was  a  pyramidal  horn  and  a  point 
source  was  taken  at  the  throat  of  the  horn  or  at  the  focal  point, 

F  (Fig.C-1 }. 

A  spherical  wave  travels  through  the  horn  and  is  incident 
on  the  excited  (curved)  side  of  the  convex  lens  {Fig.C-2},  and  is 
transformed  into  a  plane  wave. 

The  design  of  the  lens  antennas  having  plane  apertures 
(as  in  this  case)  largely  consists  of  calculating  the  profile  of  the 
excited  side  of  the  lens,  permitting  a  plane  wave  to  be  obtained  for 

given  values  of  the  refractive  index  of  the  medium,  n,  and  the  focal 

length,  f. 
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PYRAMIDAL  HORN 


(Figure  C-l) 

(all  dimensions  in  cm) 
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C.1.2  Thickness  and  Radius  of  Curvature  of  the  Plano-Convex  Lens 

The  plane  wave  in  the  lens  signifies,  firstly,  that  any  ray 
emerging  from  the  source,  incident  on  the  surface  of  the  lens, 
propagates  further  parallel  to  the  x-axis  {Fig.C-2}.  And  secondly, 
that  the  time  of  passage  of  the  ray  over  such  a  broken  path  from  the 
source  to  the  aperture  is  the  same  for  all  rays.  For  a  refractive 
index  of  more  than  unity,  a  plane  wave  characteristic  is  given  by 

1  +  2L_  =  Af  +  x)2  +  y2  /c_2 

C  V,  c  v 


and 


f  =  focal  length,  cm 

x,  y  are  the  co-ordinate  axes  {Fig.C-2} 


or 

f  +  =  /(f  +  x)*  +  y1  (C-3) 

Substitute  n  for  c/v^  from  Eq.(C-l)  into  Eq.(C-3)  to  get 
f  +  nx  =  /('f'  +  x)2  +~yT 


or 

x2(n2  -  1)  +  2fx(n  -  1 )  -  y2  =  0  (C-4) 

From  Eq . (C-4)  it  is  evident  that  for  n  greater  than  unity 
the  lens  profile  has  a  hyperbolic  form.  The  source  is  located  at  the 
focus  of  the  hyperbola  furthest  from  the  surface  of  the  lens. 

As  Eq. (C-4)  is  valid  for  any  point  on  the  curve  then  (x,  y) 
can  be  replaced  by  (d ' ,  D'/2)  the  co-ordinates  of  one  of  the 
extremeties  of  the  lens.  By  doing  so  and  solving  Eq.(C-4)  for  d  , 
the  following  equation  results. 


. 


'  - 


(D 1 /2) 2 

(n2  -  1) 
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given 


d' 


n  (for  polystyrene)  =  1.60 

f  (focal  length)  =  7.71  cm 

D 1  (aperture)  =  3.46  cm 

the  thickness  of  the  lens  d'  is  calculated  as  0.3076  cm. 

Radius  of  curvature  for  a  plano-convex  lens  can  be 
determined  by  the  relationship, 

j  =  (n  -  l)(l/Yi) 


(C-5) 


(C-6) 


i  =  radius  of  curvature. 

For  n  =  1.6  and  f  =  7.71  cm,  the  radius  of  curvature  is  calculated 
as  4.626  cm. 


C. 2  SAMPLE  CALCULATION  OF  POWER  EMERGING  FROM  THE  CORE  HOLDER 
AT  DIFFERENT  WATER  SATURATIONS 
Re-writing  Eq .III -7  as 


Ix  -4tt  f  h  k<fSw 

log10  7^  =  2.303c 


( 1 1 1-7) 


Given  that 

k  =  2.9  (value  given  in  Ch . Ill) 

<j>  =  0.37  (say) 

f  =  27 ( 1 09 )  Hz 

h  =  1.27  cm  (thickness  of  core) 

=  3(10^°)  cm/sec 


c 


' 


- 
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and  substituting  in  Eq .(II 1-7 ) 

I  4(tt)(27)(109)(1.27)(2.9)(0.37)S 

log  *  = - „ - « 

U  o  2.303(3) (1 0 1 U) 

or 

I 

lo9lOT-  =  6‘7  Sw  (C'7) 

0 

Taking  I  as  3  mW,  the  power  incident  on  the  core  holder, 
and  assuming  that  all  the  attenuation  occurs  in  the  water,  then 
Eq.(C-7)  becomes, 

^loL  =  ^ °®1 0!o  -  6-7  Sw 

or 

log10Ix  =  0.48  -  6.7  Sw  (C-8) 

In  this  simplified  calculation  for  the  power  measured, 

I  ,  the  following  values  result  when  solving  for  measured  power  at 
X 

different  water  saturations. 


TABLE  C-l 


WATER  SATURATIONS  AND  THEORETICAL 
MICROWAVE  POWER  RECEIVED 


w 

(fraction) 

0.0 

0.2 

0.4 

0.6 

0.8 

1.0 


Ix(or  R) 
(mW) 


3.00 

0.14 

0.01  4 
2.88  (10  Z) 
1.32  (10’H) 
6.03  (10_/) 


This  calculation  was  undertaken  to  have  an  idea  of  the  powers  that 
are  to  be  measured.  The  calibration  for  the  experimental  runs  was 


done  using  the  values  measured. 


. 


■ 


APPENDIX  D 

PLOTS  OF  TOTAL  PRESSURE  DROP,  PORE  VOLUMES  PRODUCED 
AND  EFFECTIVE  VISCOSITIES  VERSUS  PORE  VOLUMES  INJECTED 

DATA  IN  TABULAR  FORM  FOR  THESE  PLOTS 
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PRESSURE  (kPa) 

25.0  50.0  75.0  100.0  125 
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FIGURE  D-l  TOTAL  PRESSURE  DROP  HISTORY 
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PRESSURE  ( kPa) 

80  100  120  140  160  180  200 
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FIGURE  D-2  TOTAL  PRESSURE  DROP  HISTORY 
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PRESSURE  (kPa) 
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FIGURE  D-3  TOTAL  PRESSURE  DROP  HISTORY 
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PRESSURE  ( kPa) 
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FIGURE  D-4  TOTAL  PRESSURE  DROP  HISTORY 
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FIGURE  D-5  EFFECTIVE  VISCOSITY  HISTORY 


. 


EFFECTIVE  VISCOSITY  (mPa-sec) 

45  60  75  90  105  120  135 


FIGURE  D-6 


EFFECTIVE  VISCOSITY  HISTORY 


■ 


EFFECTIVE  VISCOSITY  (mPa-sec) 


152 


o 


« 


FIGURE  D-7  EFFECTIVE  VISCOSITY  HISTORY 
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FIGURE  D-8  EFFECTIVE  VISCOSITY  HISTORY 
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FIGURE  D-9  OIL  PRODUCTION  HISTORY 
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FIGURE  D-12  OIL  PRODUCTION  HISTORY 


A- 


TABLE  D-l 


OIL  PRODUCTION  AND  PRESSURE  HISTORY 
RUN  NUMBER  8 


T-j =  qt/ 4>AL 

Np  (cc) 

AP  (kPa) 

0.034 

6.5 

68.95 

0.068 

15.0* 

71.71 

0.109 

25.5 

79.29 

0.156 

38.0 

63.43 

0.190 

48.0 

60.68 

0.218 

55.0* 

58.26 

0.272 

68.5 

53.78 

0.340 

84.5* 

46.20 

0.408 

104.0 

37.58 

0,442 

112.0* 

35.16  (B.T. ) 

0.476 

115.5 

32.41 

0.544 

123.5 

28.96 

0.612 

126.5 

26.89 

0.694 

25.17 

1.026 

21.24 

1.109 

20.69 

* 


Saturation  profile  available 


. 


. 


TABLE  D-2 


OIL  PRODUCTION  AND  PRESSURE  HISTORY 
RUN  NUMBER  10 


x-j=  qt/<J>AL 

Np  (cc) 

AP  (kPa) 

0.021 

3.0 

33.10 

0.028 

4.5 

37.23 

0.035 

6.5* 

39.30 

0.069 

15.0* 

41.37 

0.104 

24.0* 

40.00 

0.139 

33.25 

37.58 

0.174 

42.0* 

34.48 

0.208 

50.0 

31.72 

0.243 

58.5* 

28.27 

0.278 

67.0 

24.82 

0.312 

76.5* 

21.93 

0.328 

81.0 

19.31 

0.340 

83.0* 

19.17  (B.T. ) 

0.376 

87.0 

17.93 

0.417 

91 .0 

17.58 

0.556 

103.0 

14.82 

0.694 

111.0 

12.07 

0.833 

119.0 

10.69 

0.972 

125.0 

10.00 

1.111 

9.52 

*  Saturation  profile  available 


f' 


■ 


TABLE  D-3 


OIL  PRODUCTION  AND  PRESSURE  HISTORY 
RUN  NUMBER  11 


Tj=  qt/4)AL 

Np  (cc) 

aP  (kPa) 

0.021 

3.0 

67.57 

0.048 

9.5 

95.15 

0.083 

18.0* 

96.53 

0.166 

38.5 

85.50 

0.207 

48.5* 

79.29 

0.249 

59.5 

73.09 

0.269 

65.0* 

69.98 

0.332 

71.5 

58.06 

0.363 

76.0* 

55.16 

0.390 

84.5* 

49.30  (B.T. ) 

0.415 

88.0 

45.51 

0.487 

94.0 

41.03 

0.622 

103.5 

34.48 

0.684 

107.5 

31.72 

0.767 

111.5 

29.30 

0.850 

114.0 

27.58 

0.933 

119.5 

26.20 

1.036 

122.0 

24.82 

*  Saturation  profile  available 


■ 


£ 


TABLE  D-4 


OIL  PRODUCTION  AND  PRESSURE  HISTORY 
RUN  NUMBER  12 


Ti=  qt/<f>AL 

Np  (cc) 

AP  (kPa) 

0.016 

8.0 

56.54 

0.029 

12.0* 

56.88 

0.045 

14.0 

53.09 

0.054 

16.0* 

50.68 

0.063 

18.5 

48.13 

0.072 

21.0* 

45.58 

0.082 

22.5 

43.78 

0.091 

25.0* 

41.23 

0.100 

27.5 

39.65 

0.109 

30.25* 

37.23 

0.118 

33.0 

34.82 

0.127 

35.5* 

32.96 

0.136 

36.0 

29.65 

0.145 

38.0* 

26.55 

0.152 

38.75* 

24.82  (B.T. ) 

0.163 

40.0 

23.10 

0.261 

50.9 

14.13 

0.303 

54.75 

13.10 

0.342 

57.9 

12.07 

0.386 

59.4 

11.38 

0.469 

63.4 

9.86 

0.558 

66.9 

8.62 

0.650 

70.4 

7.93 

0.688 

70.9 

7.58 

*  Saturation  profile  available 


TABLE  D-5 


OIL  PRODUCTION  AND  PRESSURE  HISTORY 
RUN  NUMBER  13 


T.j =  qt/<J>AL 

Np  (cc) 

AP  (kPa) 

0.007 

2.0 

98.60 

0.017 

4.5 

105.50 

0.034 

9.0* 

105.15 

0.051 

13.5 

99.29 

0.068 

18.0* 

92.05 

0.086 

22.5 

83.43 

0.103 

26.5* 

77.22 

0.120 

31.5 

72.40 

0.137 

35.5* 

68.95 

0.154 

40.0 

65.50 

0.171 

45.0* 

60.68 

0.188 

50.0 

55.16 

0.205 

54.5* 

49.99 

0.222 

59.5 

44.82 

0.237 

62.5* 

39.85  (B.T. ) 

0.257 

65.5 

35.85 

0.332 

73.5 

28.61 

0.397 

78.0 

24.48 

0.479 

83.0 

19.58 

0.554 

87.0 

16.89 

0.647 

90.5 

14.82 

0.708 

93.0 

13.93 

0.777 

95.5 

12.76 

* 


Saturation  profile  available 


:  p.  1 


■ 

■ 

TABLE  D-6 


OIL  PRODUCTION  AND  PRESSURE  HISTORY 
RUN  NUMBER  14 


x.j=  qt/<j)AL 

Np  (cc) 

AP  (kPa) 

O.Oll 

3.5 

141.35 

0.028 

8.0* 

142.73 

0.056 

15.0* 

129.97 

0.083 

22.5* 

105.84 

0.111 

29.5* 

97.22 

0.139 

37.5* 

82.05 

0.158 

42.0 

73.78 

0.172 

46.0 

68.61 

0.195 

51.0* 

63.43  (B.T. ) 

0.222 

57.0 

58.61 

0.278 

67.0 

50.33 

0.334 

75.0 

37.23 

0.417 

83.0 

27.24 

0.528 

87.5 

21.72 

0.612 

92.0 

18.96 

0.695 

95.0 

16.89 

0.778 

98.0 

15.51 

0.862 

100.5 

14.48 

0.945 

103.0 

13.45 

* 


Saturation  profile  available 


TABLE  D-7 


OIL  PRODUCTION  AND  PRESSURE  HISTORY 
RUN  NUMBER  15 


x.j  =  qt/4>AL 

Np  (cc) 

AP  (kPa) 

0.034 

9.0 

193.06 

0.069 

19.0 

174.44 

0.103 

28.5* 

155.83 

0.138 

37.5* 

136.52 

0.172 

46.5* 

116.53 

0.207 

54.0* 

97.91  (B.T. ) 

0.275 

67.0 

68.95 

0.344 

75.0 

53.09 

0.413 

81 .0 

43.09 

0.489 

86.0 

37.23 

0.558 

90.0 

33.79 

0.661 

94.0 

29.65 

0.792 

99.0 

25.51 

0.895 

102.5 

23.44 

* 


Saturation  profile  available 


TABLE  D-8 


OIL  PRODUCTION  AND  PRESSURE  HISTORY 
RUN  NUMBER  16 


t. =  qt/fAL 

Np  (cc) 

aP  (kPa) 

0.007 

1.5 

26.20 

0.036 

8.5 

33.10 

0.072 

16.5* 

32.75 

0.109 

25.0 

32.06 

0.145 

32.5* 

31.72 

0.181 

41.0 

31.37 

0.217 

50.5* 

30.34 

0.254 

59.0 

28.96 

0.290 

67.0* 

26.20 

0.326 

76.0 

23.44 

0.362 

83.5* 

21 .03 

0.399 

93.5 

18.62 

0.435 

102.5* 

15.86 

0.460 

106.5* 

13.79  (B.T. ) 

0.507 

112.0 

13.10 

0.638 

120.0 

10.34 

0.725 

124.75 

9.65 

0.798 

128.25 

8.96 

0.870 

131.25 

8.62 

0.942 

134.25 

8.27 

1.014 

137.25 

7.93 

*  Saturation  profile  available 


* 

. 

TABLE  D-9 


OIL  PRODUCTION  AND  PRESSURE  HISTORY 
RUN  NUMBER  18 


T.j  =  qt/cfAL 

Np  (cc) 

AP  (kPa) 

0.007 

0.6 

6.21 

0.036 

5.5 

37.58 

0.072 

14.0* 

39.30 

0.109 

22.5 

38.61 

0.145 

30.5* 

37.23 

0.181 

39.5 

35.85 

0.221 

48.5* 

33.79 

0.254 

56.5 

31.72 

0.290 

65.0* 

28.96 

0.326 

74.0 

26.20 

0.362 

83.0* 

22.75 

0.399 

91.5 

20.00 

0.435 

99.5* 

17.24 

0.449 

102.0 

16.55  (B.T. ) 

0.471 

107.5 

15.17 

0.507 

112.5 

13.79 

0.587 

118.0 

12.41 

0.652 

123.0 

11.38 

0.725 

127.0 

10.34 

0.797 

130.2 

9.65 

0.870 

132.45 

9.31 

0.942 

135.45 

8.62 

1.014 

138.45 

8.27 

*  Saturation  profile  available 


. 


TABLE  D-l 0 


OIL  PRODUCTION  AND  PRESSURE  HISTORY 
RUN  NUMBER  19 


Tj=  q  t/<f  AL 

Np  (cc) 

AP  (kPa) 

0.014 

2.5 

28.27 

0.036 

7.5 

33.79 

0.072 

16.5* 

32.41 

0.109 

25.0 

31.03 

0.145 

33.0* 

29.65 

0.181 

43.0 

27.24 

0.217 

51.5* 

24.82 

0.254 

59.0 

22.75 

0.290 

68.5* 

20.34 

0.326 

78.0 

17.93 

0.362 

86.5* 

15.51 

0.378 

90.5 

14.48  (B.T. ) 

0.399 

94.5 

13.79 

0.435 

98.5 

12.41 

0.507 

105.5 

10.86 

0.580 

110.75 

9.65 

0.652 

115.5 

8.96 

0.725 

119.5 

8.27 

0.797 

123.5 

7.58 

0.873 

127.0 

7.24 

0.942 

130.25 

6.90 

1.014 

133.75 

6.76 

*  Saturation  profile  available 


TABLE  D-l 1 


OIL  PRODUCTION  AND  PRESSURE  HISTORY 
RUN  NUMBER  20 


ij  =  qt/(j)AL 

Np  (cc) 

aP  (kPa) 

0.072 

12.5* 

67.57 

0.145 

29.5* 

64.81 

0.217 

47.5* 

58.61 

0.304 

65.5* 

48.27 

0.362 

83.0* 

40.68 

0.394 

91.0 

35.85  (B.T. ) 

0.435 

95.5 

32.41 

0.507 

105.5 

26.89 

0.580 

113.5 

23.44 

0.652 

119.0 

21.37 

0.725 

123.0 

20.00 

0.797 

127.0 

18.96 

0.870 

129.5 

17.93 

0.942 

132.0 

16.89 

1.014 

134.5 

16.41 

* 


Saturation  profile  available 


•  • 

- 
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TABLE  D-12 


OIL  PRODUCTION  AND  PRESSURE  HISTORY 
RUN  NUMBER  21 


Tj  =  qt/(j)AL 

Np  (cc) 

AP  (kPa) 

0.109 

23.0* 

108.25 

0.217 

48.0* 

85.50 

0.326 

76.0* 

65.50 

0.397 

94.0 

53.10  (B.T. ) 

0.435 

100.0 

48.27 

0.543 

112.0 

37.23 

0.652 

118.5 

32.41 

0.761 

124.5 

29.65 

0.870 

129.5 

26.89 

0.978 

134.0 

25.17 

1.087 

137.5 

23.44 

* 


Saturation  profile  available 


TABLE  D-l 3 


OIL  PRODUCTION  AND  PRESSURE  HISTORY 
RUN  NUMBER  9 


-e- 

+-> 

C r 

II 

H- 

Np  (cc) 

aP  (kPa) 

0.014 

1.0 

4.14 

0.070 

15.0 

11.38 

0.140 

32.0* 

13.45 

0.210 

49.0* 

15.51 

0.245 

58.0 

16.55 

0.320 

76.5 

18.75 

0.356 

85.75* 

20.00 

0.420 

102.0 

21.86 

0.456 

110.0* 

23.10 

0.525 

127.5* 

25.86 

0.600 

145.0 

27.58 

0.630 

153.5 

29.10 

0.644 

154.5 

30.34  (B.T. ) 

* 


Saturation  profile  available 


TABLE  D-14 


OIL  PRODUCTION  AND  PRESSURE  HISTORY 
RUN  NUMBER  23 


x-j  -  qt/(J)AL 

Np  (cc) 

AP  (kPa) 

0.036 

6.5 

8.96 

0.080 

16.5 

10.34 

0.109 

23.5 

10.69 

0.146 

31.5* 

10.89 

0.189 

41.0 

12.07 

0.218 

48.5* 

12.41 

0.255 

57.5 

14.13 

0.290 

65.5* 

15.17 

0.328 

75.5 

17.24 

0.364 

83.5* 

18.62 

0.400 

92.5 

20.69 

0.437 

109.5* 

22.06 

0.473 

119.0 

23.79 

0.509 

128.0* 

25.86 

0.546 

136.5 

27.58 

0.582 

145.5* 

31.37  (B.T. ) 

* 


Saturation  profile  available 


• 

TABLE  D-l 5 


OIL  PRODUCTION  AND  PRESSURE  HISTORY 
RUN  NUMBER  24 


T-j  =  qt/<j)AL 

Np  (cc) 

aP  (kPa) 

0.036 

6.5 

5.17 

0.072 

16.0* 

5.38 

0.109 

24.5 

5.52 

0.145 

32.5* 

5.52 

0.181 

42.0 

5.64 

0.217 

51.0* 

5.79 

0.254 

58.5 

6.90 

0.290 

67.0* 

7.03 

0.308 

72.0 

7.24 

0.326 

76.0 

8.27 

0.362 

84.5* 

8.62 

0.399 

93.0 

9.65 

0.435 

o 

• 

o 

* 

10.34 

0.471 

110.5 

11.03 

0.507 

119.5* 

11.72 

0.543 

128.0 

12.76 

0.587 

136.5 

16.55  (B.T. ) 

* 


Saturation  profile  available 


' 

APPENDIX  E 

EULERIAN  SOLUTION  PROCEDURE  FOR  IMMISCIBLE 
DISPLACEMENT  EQUATION 
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E.l  EULERIAN  SOLUTION  PROCEDURE  FOR  THE  IMMISCIBLE  DISPLACEMENT 


EQUATION 

Consider  Eq.(II-7),  and  differente  both  sides  with 
respect  to  c  to  have 

3f. 

v  _  c 

dS*  W 


9? 


w  .  -  JL  {ncC(S*)  If } 


(E-l) 


From  Eq.(E-l),  substitute  the  value  of  af/ac  in 

w 


9? 


{Ncc<s*>  If}  - 


dG(S*)  as* 


Eq. ( I I -6 ) ,  to  get 

as* 

ax*  aU  Vw  '  a?  /  dS*  ac 
One  way  of  handling  the  second  order  derivative  in 
Eq.(E-2),  with  the  non-1 ineari ty ,  is  to  use  the  differencing  scheme 
proposed  by  Blair  and  Weinaug^  [1969,  p . 41 8] , 


(E-2) 


_a_ 

9? 


{n  c(s*)  —  \  =  — 
l  c  v  1  ac  f  ac 


s*  -  S*)  (S*  -  s*  i 

w  <  Jit1. z)  n  n_1 


AC 


AC 


( E-3) 


where 


w  =  f{c(S*)  +  C(s*+1)} 


(E-4) 


and 


N 


Z  = 


_c 

2 


{C(Sn-l  >  +  C(Sn)| 


(E-5) 


A1  so 


AC  =  step  size  for  the  space  variable 
n  =  counter  for  the  spatial  nodes. 


Substituting  for  ^-|ncC(S*)  ||—j  from  Eq . (E-3) 
one  obtains, 

as*  .  dG ( S* )  as*  1  Li  S?i+1  ■  Sn  } 
a?F  -  alT  AsL  j~~A5  j  L 


into  Eq. (E-2) 


(S*  -  S*  , 
)  n  n-1 


(E-6) 


. 


*  £ 
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Considering  i  as  the  counter  for  points  in  time  and 


di scretizing 

the  remaining 

derivatives  in 

Eq.(E-6)  to  have 

$*i 

-  S*i-1 
n 

dG(S*i) 

1  ( 

Sn+1  -  Sn-1  } 

Ax* 

dS* 

2AC  \ 

♦  w  j 

(AS)2  1 

(  c*i 
{  n+1 

‘  S^} 

-  Z  \ 
(AS)2  ' 

hs1  - 

S^l} 

(E-7) 

V 

Simplifying  Eq.(E-7)  to  obtain, 

-  (A+E)S*i-, 

+  (l+B+EjS*1 

+  (A-B)S; 

*i 

n+1 

c*  i  - 1 
n 

( E-8 ) 

where 

A 

dG(S*i) 

“  m  m  • 

At* 

( E-9 ) 

dS* 

2AC 

B 

=  WAT* 

(as)2 

(E-10) 

and 

E 

■7  At* 

(as)2 

(E-n) 

Mote  that  Eq.(E-8)  is  fully  implicit  in  the  time  derivative. 
Saturations  at  (i-l)th  time  step  are  known  and  are  to  be  evaluated 
for  the  i-th  time  step. 

E.1.1  Initial  and  Boundary  Conditions 

The  conditions  used  in  the  numerical  procedure  were, 

S*(s,  0+)  =  F(s)  0  <  S  s  1  (E- 12) 

S*(l,  T*)  =  0  0  <  x*  <  ( E- 1 3 ) 

and  to  get  the  numerical  procedure  started  the  inlet  boundary 
condition  for  the  first  time  step  is 


. 


. 


. 
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S*(0,  At*)  =  SI  (E-14) 

where  SI  is  assigned  an  arbitrary  value  between  zero  and  one.  The 
Eq.(E-8)  is  solved  with  Eqs . ( E- 1 2 )  to  (E-14)  for  the  first  time  step 
and  material  balance  is  checked  by  numerically  integrating  for  the 
area  under  the  first  profile.  The  area  must  be  equal  to  At*,  the 
magnitude  of  the  time  step,  if  the  choice  of  SI  is  correct  for  that 
particular  data  set. 

F(s)  is  also  arbitrary,  and  in  this  numerical  procedure  a 
convenient  choice  was, 

F(c)  =  -  —  5  +  SI  C  i  AC 

or 

FU)  =  0  5  >  A? 

It  is  evident  that  F(?)  is  different  for  every  choice  of  SI. 
In  actuality,  this  function  is  only  used  in  the  first  iterate  and  then 
the  current  nodal  values  are  used  till  convergence  is  achieved. 

Once  the  profile  after  the  first  time  step  is  obtained,  the 

saturation  at  the  inlet  will  build-up  for  the  subsequent  time  steps. 

To  enable  the  saturation  to  build-up  Eq . ( E- 1 6 )  was  used.  Once  again 

differentiating  both  sides  of  Eq . ( 1 1-7 )  with  respect  to  £  and 

substituting  the  resulting  3f.,/3?  in  Eq .  ( 1 1 -6 )  to  obtain, 

w 

3S1  -  N  <  3C (S*)  3S^  +  jfS*  _  3G(,S*),  , E, 6) 

3t*  'M  3C  3C  'f  3C  U  ' 


(E-15) 


The  finite  difference  analog  of  Eq.(E-16)  for  node  1  can  be  written  as. 
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S*i+1 

1 


NcAt* 

(A?)2 


|c(S*' )  -  C ( S* 1 ) J-  (S*1  -  S*i) 


N  At* 


{s?1  - 


2S*1  +  S*1 


v) 


GfS*1)  -  GtS*1) 


AC 


+  s*1 


(E-17) 


Once  the  inlet  end  saturation  attains  a  value  of  unity  the 
third  boundary  condition  becomes, 

S*(0,  T*)  =  1 

Eq.(E-8)  was  written  for  every  node  and  the  resulting 
coefficient  matrix  was  tri -diagonal ,  facil ititating  solution  by 
Thomas's  procedure.  For  every  time  step  S*1' ,  C(S*i)  and  dG(S*i)/dS* 
were  iterated  upon  till  convergence  was  achieved.  The  time  step  size. 
Ax*,  was  taken  as  0.005  while  the  spatial  axis  was  divided  into  one 
hundred  parts  to  give  a  step  size,  AC,  of  0.01. 


